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Abstract 

Contrary to the finite dimensional case, Weyl and Wick quantizations are no more asymptotically 
equivalent in the infinite dimensional bosonic second quantization. Moreover neither the Weyl calculus 
defined for cylindrical symbols nor the Wick calculus defined for polynomials are preserved by the action 
of a nonlinear flow. Nevertheless taking advantage carefully of the information brought by these two 
calculuses in the mean field asymptotics, the propagation of Wigner measures for general states can be 
proved, extending to the infinite dimensional case a standard result of semiclassical analysis. 

2000 Mathematics subject classification: 81S30, 81S05, 81T10, 35Q55 

1 Introduction 

Our main result is briefly presented in this introduction. Accurate definitions will be found in Section|2] 
Let Jf? = Y s (3f ) be the bosonic Fock space constructed over the complex separable Hilbert-space iF, 
T s {3f) = (£>™ =0 \/ n 3f where V" ^ is the symmetric n-th hilbertian tensor power of 3? . Consider the 
Hamiltonian 

ff e = dT(A) + (£ (z® J ,Qjz® } )) Wick 

7=2 

defined for the self-adjoint operator (A, &{A)) on J° and Qj = Qj* G -^(V ; It is the Wick quantized 
version of the classical Hamiltonian 

h( z ,z) = (z,Az) + j^ (z^,QjZ® j ) , z e 9{A) c 2. 
When 3f =L 2 (R d ) , the operator H £ is formally written 

H £ = A(x,y)a*(x)a(y) dxdy+Y / Qj(xi,...,Xj,yi,...,yj)a*(xi)..M*(xj)a(yi)...a(yj)dxdy, 

with the £-dependent canonical commutation relations [a(x) , a* (y)] — e8(x — y). Here A(.,.) and Qj(.,.) 
denote the kernels of the operators A and Qj. The mean field asymptotics is concerned with the limit as 
e 0, where ^ = N £ represents a large number of particles and where e enters in the CCR-relations by 

Vf,g€2r, [a(f),a*(g)] = e(f,g)I. 
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The number operator is N = dF(L r ), with Nz®" = Enz m . For a normal state p £ £ Jz? 1 (V^ ^) C 1 ( ^) 
with = L 2 (R d ), a standard tool considered in the mean field limit is the BBGKY hierarchy of reduced 
density matrices: 

$' ) (x,y)=[ Pe (x,X,y,X)dX, peN, 

J R 2d(N £ -p) 

and such a definition will be extended to general 2f and normal states p £ £ Jz? 1 (J$?) fulfilling the condition 
Tr [p e N*] < +°c for all A: € N . 

For a cylindrical function, b(z) = b(pz) for some finite rank projection p and b belonging to the Schwartz 
class S^{p3f), the Weyl quantization can be given by 

b Weyl = j ^ (z) w ( Lp ^ ^ 

where W(V2?rz) = e' n ^ z ' +a and where L p and & are respectively the Lebesgue measure on p2£ and 
the (e-independent) Fourier-transform on ,5^(p3f ). Associated with a family (pe) ee (0,e)' Wigner measures 
can be defined by 



lim Tr 



after extracting subsequences under the sole uniform estimate Tr [p £ N 5 ] < Cg for some 5 > . 

The problem of the mean field dynamics questions whether the asymptotic quantities as e — > associ- 
ated with 

p e (t)=e- i ° H *p e e-^ H *, t£R 

are transported by the flow F, generated by the classical Hamiltonian h(z,z) and given, after writing zt = 

V,- S (z s ),by 

id,z, = m(z t ,zr) =Az t + £ j(z? ] - X , QjZ? ] ) . (1) 

7=2 

The finite dimensional case enters in the standard framework of semiclassical analysis and has been studied 
extensively in the 80's and 90's by various authors and with various methods ( SUE] Ell E2 Ifl8llll43llll24ll 
and references therein). 

It was first considered by Hepp in ll36l and extended by Ginibre and Velo in [30] [31] by the squeezed 
coherent states method well-known as the Hepp method (see also [49||6|). More recently the question 
of the mean field dynamics has been tackled with the so-called BBGKY-hierarchy approach inspired by 
the BBGKY-method of classical kinetic theory (see J52) (Bl E2 C31 EUffl EllESland also the related 
works ll37llfT7l ). In |25||26||27| a specific use of the structure of the Wick calculus in the bosonic Fock 
space was used to make work truncated Dyson expansions for the mean field dynamics of specific states. 
The aim of our work started in Q was to restore the phase-space geometric nature of the problem in the 
spirit of Ull [33 ]| 38 ]| 39 1 and to extend as much as possible to the infinite dimensional case, the methods 
well understood for the semiclassical finite dimensional problem. In this first article, we explained the 
construction of Wigner measures, analyzed accurately the gap of information carried by Weyl observables 
and Wick observables and use these Wigner (or semiclassical) measures to reformulate known propagation 
results. In [8], we reconsidered the truncated Dyson expansion method of [25 1[ 26] 1 27] in order to prove the 
propagation of Wigner measures for some specific families of states. We are now able to state the following 
general result (still with a regular interaction term contrary to many other works cited above). 

Theorem 1.1 Let (ps) e e(o,e) be a family of normal states on Ji? with a single Wigner measure /j.q and such 
that 

VaeN, limTr[p e N a ]= / |z| 2a d^{z) < +°°. (2) 

Then for all t £ M, the family (pe(f) = ^' lHe Pe^' lHe )Ee(0.e) has a unique Wigner measure jJ, t = (F^/io, 
which is the initial measure jUo pushed forward by the flow associated with ([T). 
Moreover the convergence 



lim Tr 

£^0 



p £ (t)b Wlck = j bo¥ t (z)d^(z) 
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holds for any b € @> alg {2e) = ®J^9 PA {&) ■ 
Finally, the convergence of the reduced density matrices 

holds in the Jzf 1 (V P T)-normfor allpEN. 

Comments: The existence of Wigner measures as Borel probability measures requires a uniform estimate 
Tr [peN 5 ] < C§ for some 8 > 0, but such an assumption would be redundant with the existence of bounded 
limits stated in (0. 

The uniqueness of the Wigner measure is not really a strong assumption since it suffices to replace the 
whole family (pe)ee(0,e) by a suitable extracted sequence (pe k )keN, timk-^°° e k = 0, in order to fulfill this 
requirement. Such a reduction argument after extraction will often be used. 

The fact that the quantities Tr[p £ N a ] are uniformly bounded w.r.t £ 6 (0,e) is also very natural within the 
mean field framework and satisfied by all known physical examples. 

Actually the strong assumption which is not satisfied in all cases is that the limit in (O equals \z\ 2a djM)- 
This condition prevents from the phenomenon of "infinite dimensional defect of compactness" identified 
in ||7'| and which was shown to appear in the physical example of the Bose-Einstein free gas (the non 
condensated phase is responsible for a discrepancy between the left- and right-hand sides of (O). The 
analysis of this phenomenon is improved in Section [2] 

Finally our proof no more uses truncated Dyson expansions of the quantum flow and relies only on the 
good properties of the classical flow, after exploiting all the a priori information given by the Weyl and 
Wick calculus. 

Outline: The Section [2] introduces the various objects used for our analysis, Wick and Weyl calculuses, 
Wigner measures, reduced density matrices. The conditions presented in |8| are reduced to the simple 
equivalent form (fJJ in Subsection l2.7l After this the Subsection l2.8l is devoted to the notion of states local- 
ized in a ball. 

The dynamics is studied in Section[3] First a simple condition is proved to ensure, via some equicontinuity 
argument, the possibility of a common extraction process (s^keN f° r a ^ times t € R . Then the propagation 
of Wigner measures is proved for states localized in a ball. Then the truncation is removed and all the argu- 
ments are gathered for the proof of Theorem l 1 . 1 l in Subsection l3.4l Finally, additional simple consequences 
are listed in Subsection l3.5l 

Examples are presented in Section@] It is recalled that the regular interactions are physically relevant within 
the modelling of the rapidly rotating Bose condensates in the Lowest Landau Level approximation. Details 
are given about the propagation of non trivial Wigner measures supported on a torus, which shows the ad- 
vantage of this formulation compared to the BBGKY hierarchy method. Finally, the propagation of Wigner 
measures provides a nice formulation of the Hartree-von Neumann limit. 



2 Information carried by Wigner measures 

After introducing the symmetric Fock space with £-dependent CCR's and recalling some properties of 
the Wick quantization, the connection between infinite dimensional Wigner measures and the BBGKY 
presentation of the many body problem is explicitly specified. This section ends with the notion of states 
localized in a ball, which will be usefull in the proof of Theorem ll.il 

2.1 Fock space 

Consider a separable Hilbert space 2f endowed with a scalar product (., .) which is anti-linear in the left 
argument and linear in the right one and with the associated norm \z\ = sj (z,z). Let a = Im(.,.) and 
5 = Re(.,.) respectively denote the canonical symplectic form and the real scalar product over 3f. The 
symmetric Fock space on 2f is the Hilbert space 

^=0Y»f = r,(f), 

n=0 
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where V" 3? is the n-fold symmetric tensor product. Almost all the direct sums and tensor products are 
completed within the Hilbert framework. This is omitted in the notation. On the contrary, a specific alg 
superscript will be used for the algebraic direct sums or tensor products. 

For any n G N, the orthogonal projection of ®" H£ onto the closed subspace V" will be denoted by 
y n . For any , £>, ...,£,) G the vector |i V | 2 V • • • V £„ G V" # will be 

^iV§ 2 V---V^ = t 5^(§i®fe---®^) = -^ E &r(i)®6c(2)-"®&(»)' (3) 

where 6„ is the symmetric group of degree n. The family of vectors (|i V • • • V ^ n )^e2" is a tota l family of 
V" ^ and thanks to the polarization identity 

SiV&v-v^-L £ ei .-. e „(f e^-r, (4) 

£,=±1 .7=1 

the same property holds for (4 8 ") n eN • 

For two operators A k : V'* 3f -> V 7 * ^,£=1,2, the notation A i V A 2 stands for 

Ai V^2 = ^ 1+/2 o (Ai ®A 2 )oy h+h e jgf(V h+h ar,y h+h ar) . 

Any ze Jis identified with the operator |z) : V° ^ = C 9 A Xz G = V 1 & while (z| denotes the 
linear form f 9^(z^)eC. The creation and annihilation operators a* (| ) and a (| ), parametrized by 
£ > 0, are then defined by: 

«(4)|V";r = Ven {Zl^Iyn-ig. 

a*(^)|V"^ = A/ £ ( n + 1) ° ( I4> ® V*) = vM" + 1)4 VViT 

and satisfy the canonical commutation relations (CCR): 

[fl(4i),fl(&)] = [fl*(4i),fl*(&)] =0, [«(&),«*(&)] =£<«!,&>/. (5) 
We also consider the canonical quantization of the real variables <E>(4) = («*(<!;) + a(4)) and 11(4) = 
<£>(*'<[;) = -^( a (4) ~ a *(4))- They are self-adjoint operators on and satisfy the identities: 

[*(&),*(&)] = fea(4i,&)7, [*(5i),n(&)] = i£S(4i,&)/. 
The representation of the Weyl commutation relations in the Fock space 

w(4i)w(4 2 ) = e -¥ CT( ^w(4i + 42) (6) 

= <r ie<T ^%(§ 2 M§i), 
is obtained by setting W(4) = e'*^. The number operator is also parametrized by e > 0, 

N| y« j» = £n/| yn tsg . 

It is convenient to introduce the subspace 

nGN 

of , which is a set of analytic vectors for N. 

For any contraction S G j£f (<2Q, \S\^^^ < 1, T(5) is the contraction in defined by 

T{S)\ v ,2° = S ® S 1 • • • (g> 5 1 . 

More generally r(B) can be defined by the same formula as an operator on J^,„ for any B G Jf(3f). 
Meanwhile, for any self-adjoint operator A : 2f D 3>{A) — > iF, the operator dr(A) is the self-adjoint oper- 
ator given by 

e |dT(A) =r(e , t A ) 



dT(A)| V „,ai g ^ (A) - e 

For example N = dr(7) . 



£/<g)---<g)^A 

*=1 



ft 
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2.2 Wick operators 

The Wick symbolic calculus on (homogenous) polynomials as introduced in [7 1 is recalled with its basic 
properties. 

Definition 2.1 For p,q £ N, £Pp„(3f) denotes the set of (p,q)-homogeneous polynomial functions on 3? 
which fulfill : 

b{z) = {z m M®») with b£%(S] p 2?\j q 2?)- 

The subspace of & p , q {3?) made of polynomials b such that b is a compact operator b € Jz?°° (V P ^ i V" 7 -2°) 
(resp. be^ r {\l p ^,\l q ^)) is denoted, by S^ A {St) (resp. @> r M {T)). 
On those spaces, the natural norms are 



\ b \^ P .q = l^ljf(V p zwa?) and \ b \& r p , q = \ b W{\ p i 1 < r - 

\ al S (St ( Oi 



The set ofnon homogeneous polynomials, the algebraic direct sum ® a J S ae m^'p,q(^') ( resp. 9°'%^! q{^) 



with 1 < r < °°), will be denotedby ^algi^) (resp. 0*^ (&)). 

Owing to the condition b € J£(\/ p J", \/ q for b € ^ Pjq {3f), this definition implies that any Gateaux 
differential dld k b(z) at the point zef belongs to jSf (V k 3T, V > ' with 

(cp,did*b(z)Y) = p \ j! .,, (^V9,^- t Vy) . 

[p-ky. {q-jy. 

In particular, we recover the operator b from b(z) via the relation 

5 = II e -Sf (W V 9J? ) ■ 

p\ q ] . v v 

With any "symbol" b € S? p ^(^\ a linear operator b Wtck called Wick monomial can be associated according 
to: 

b Wlck : Jiff i n -» ^/ in , 



with /? = (q\)~ l d?d q b(z) . The basic symbol-operator correspondence: 

(z,£) <— ► a* (4) V5S(£,z) <— ► <&(£) (z,Az) <— ► dT(A) 

(|,z) 4— ► fl(|) v^a^z) n(S) |z| 2 <— ► N, 

and more generally 

fn^^xn^z) j =fl*(i?i)-a*(i7x§i)-fl(§ ff ). 

We have the following properties. 

Proposition 2.2 The following identities hold true on Jiff for every b € ^ pq (^): 
(i) (b wick Y =W ick . 



(ii) (C{z)b(z)A(z)) Wk ' k = c Wick b Wkk A Wick , if A e & afi (2r), C e ^ojsC^)- 
(iiij e ^dr(^)foW^ e -i|dr(A) = (i,( e -M z )j« ; y A - ? fl se if. a dj int operator on 3f. 

A consequence of i) says that b is symmetric when g = p and = b. Moreover the definition (0) gives 

(q = p and b > 0) (& w/c * > on J^„) , (8) 
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which is false for general non negative polynomial symbolfl For an increasing net of non negative opera- 
tors (b a ) a , b a G J£(\/ p 3f) (again q = p), it also gives 

(b = sup b a in J2f(V P ^)) => fv<p G Jfyk , (9 , fc m ' c *<p) = sup(<p , ft™*?)) • (9) 

When & = L 2 (R d ,dx), the general formula for b wick with b G Pp^ST) is simply 

jFk* = / fe^,..., y 9 ,jci,...,Xn)a*(yi)...a*(yo)fl(jci)...a(jCp) dx\---dx p dy\---dy q , 

where is the Schwartz kernel of b and where a(xk) — a(8 Xk ) according to the usual convention. 

Proposition 2.3 For b G £P > p ^(£tf), the following number estimate holds 

(N)-h w,ck (N)-? <\b\& . (10) 

The relations © and © now become for b G @> p , p {3?) or fc a G £^(3?) 

(<? = /;> and £>0) => ((N)- p / 2 fc Wic/: (N)- p / 2 >0 inJSf(^f)) , (11) 
(5 = su Pa £ a in ^(V p 5T)) =>> ((^-"/^^(N)-^ 2 = su Pa (N)- p / 2 bf ck (N)- p / 2 in jSf(jg°)) .(12) 

An important property of our class of Wick polynomials is that a composition of b^' ck o b^' ck with b\ ,b 2 G 
Palgi^O is a Wick polynomial with symbol in &> alg (T). For G ^ pim (3T), b 2 G ^ P2)92 (^') ) G 
N and any fixed ot**i(z) G jSf(V*^;C) while d 2 *6 2 (z) G \J k 3T. The C-bilinear duality product 

d k b\(z).d k b 2 {z) defines a function of z G 5° simply denoted by d k b\.d k b 2 . We also use the following 
notation for multiple Poisson brackets: 

{fti,*2} (t) =d k bi.d k b 2 - d k b 2 .d k b u keN, 
{b h b 2 } = {b h b 2 } W . 



Proposition 2.4 Letb\ £ &> puqi (2?) andb 2 G ^ nm ^). 

For any k G {0, . . . ,min{pi ,^2}}. <5|^1 -<5|^2 belongs to 3 g P2 -k,q l -k(^') with the estimate 



1 z i z m. Sp2.cn — / £\ J / k)l L ' i ^Pi'ii< z ' ir P2>n 

The formulas 

(■ r 1 \ Wick 

mm{ P1 ,q 2 } k \ \ Widt 

I fi ^ 1 '^ 2 ) = (e e(5z,<?ffl> M«»)| z = 

/max{min{p 1 ,,jr2},min{p2,9l}} \ Wic* 

(/0 [C^^> E ^{W fc) , 



/10W as identities on Jtff;, 



fin- 



2.3 Cylindrical functions and Weyl quantization 

Let P denote the set of all finite rank orthogonal projections on 3? and for a given p G P let L p (dz) denote 
the Lebesgue measure on the finite dimensional subspace p^f. A function / : 3? — > C is said cylindrical if 
there exists p G P and a function g on p^f such that f(z) = g(pz), for all z G In this case we say that / 
is based on the subspace p2£ ' . We set .5"%/ (3°) to be the cylindrical Schwartz space: 

(/ G y cy ,m) (3p G P,3g G J^(/>iF), /(z) = g(pz)) . 

1 This property should not be confused with the positivity of the finite dimensional Anti-Wick quantization which associates a non 
negative operator to any non negative symbol. 
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The Fourier transform of a function / G based on the subspace p3f is defined as 

Jn3f 



and its inverse Fourier transform is 



/(*)=/ -n/Kz) e 2 * iS ^L p (dz). 



With any symbol b € y cy i(3f) based on p3f, a Weyl observable can be associated according to 

b Weyl = I jr[ fo ]( z ) W(V27tz)L p {dz). (13) 



After the tensor decompositions 



= r s (2f)=r s {p2?)®r s ((l-p)3?) due to 3T = p2r®(l-p)2f 

Vz € p3T, W(V2tzz) = W pS {V27tz)®Ir si i- P )2- 

where W p & denotes the reduced representation in F s (p2f), one sees that the Weyl quantization of cylin- 
drical observables based on piF amounts to the usual finite-dimensional Weyl quantization. Hence more 
general classes of symbols can be considered. 

For p £ F, the symbol classes defined for < V < 1 on the finite dimensional phase space p2£ , 

= ™-), (14) 

where {z) 2 p = 1 + are natural Weyl-Hormander algebras associated with the finite dimensional har- 

monic oscillator Hamiltonian, N p = (\z\ p3f ) Wick = (\z\ 2 p3f ) Weyl ~ ^^-£. They contains the polynomial 
functions on p2f. The associated class of Weyl quantized operators after tensorization with IrJ(\-p))% is 
denoted by OpS^. For a cylindrical polynomial b e £P a i g (3f) based on p3f, b{z) = b(pz), the asymptotic 
equivalence of the Weyl and Wick quantization in finite dimension says for any v € [0, 1] 

b Wick = b weyi + ffh ^ in QpS v^ (15) 
Such polynomials have finite rank kernels and make a dense set in but not in {P a i g {3F). 

2.4 Wick observables and BBGKY hierarchy 

When 2? = L 2 {R d ), mean field results are often presented or even analyzed in terms of reduced density 

matrices or more precisely in terms of a sequence {% ) p eN with Ye G Jz? 1 (V p -2°) ■ This follows the 
general BBGKY approach of the kinetic theory and the y£ correspond in the classical case to the empirical 
distributions. 

The basic example is when p £ e -S? 1 (V" « = [\] ■ For any p eN,p< n, yf ] G if 1 (V p &) is defined 
as the partially traced operator with the kernel 

ji p \x u ... : x p ;yi ...y p ) := / p e (x 1 , . . . ,x p ,X,y u . . .y p ,X) L^Ha-p) {dX) . 

With the polarization identity ©, the family {\w m ){W m \)w^^ f° rms a total set of jSf^V" 2?). Hence the 
formal identity 

£P jSp)\ I V| 2( '" P V(^i) ■ • • W(X P )WM ■ ■ -Wp) = (a(yi). ..a{y P )w m , . ..a(* P )y^ B ) 

= Tr [a* Cv! ) . . . a*(y>(*i ) . . . a(x p ) \ O ( y®" |] 

carries over to p £ e Jz? 1 (V" : 

Vpe {!,...,«} , £P J^—py ri P) (x u ■ ■ ■ ,x p ,yi . . . ,y p ) = Tr[a* ( yi ) . . .a* (y p )a(xi) . . .a(x p )p £ ] . 
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The correct meaning of this definition is 



Tr 



i p) h 



£Pn(n- l)...{n- p+l) 



Tr 



Wick 



Vb££? PiP (2f). 



Moreover after noticing that the factor e p n(n — l)...(n — p+ 1) is nothing but Tr [pedzl 2 '') 1 ^'*] when 
Tr [p e ] = 1 and p £ e Jgf 1 (V" it becomes 



Tr 



#5 



Tr[p e 



Tr[p £ (k| 2 ") WM 



■Tr 



Wick 



(16) 



with the convention that the right-hand side is when Tr [pe(|z| 2p ) W ' !<: *l = 0. The extension to general 
p £ £ Jz? 1 (Jf?) requires an assumption. Moreover it works for a general separable Hilbert space 3f. 

Proposition 2.5 Assume that p £ £ ££ 1 satisfies p £ > andN k / 2 p £ N k/2 £ Jz^ 1 (Jtf )for all ken. Then 
for any p £ N, f/ie relation ( fT6l l defines a unique element yi P) > o/J^ 1 (V p ^) . 

Proof. Suppose Tr [p £ (|z| 2p ) w, ' cA 'l > 0. Writing 



Tr 



Pefc 



Tr 



(1 + N) p/2 p £ (l +N)"/ 2 (1 +N)- p / 2 b Wick {l +N)- p / 2 



with our assumptions and the estimates (TlOb ensures that & — > Tr [p £ Z? W!C *] defines a continuous linear form 
on Jz? (V P • The positivity comes from (fTTT i and the normality of the associated state after normalization, 
which says yf ] £ if 1 {\J P 2f), is a consequence of i\2\ . □ 



We end with this discussion with a natural definition. 



Definition 2.6 When p £ £ Jz? 1 (&) .sato/es p £ > ant/ N^pgN*/ 2 € Jz? 1 (^) /or aH iteN, f/ie ret/weed 
density matrix jf\ pgN, associated with p £ is the element of Jz? 1 ( V P -2°) defined by 



Tr 



^5 



Tr[p £ 



Tr[p f (|zpP)W^ 



■Tr 



P £ £ 



(17) 



H-iYA }^ p) = in fne case w/ien Tr [p £ (|z| 2p ) w, ' c *] = 0. 



2.5 Wigner measures 

The Wigner measures are defined after the next result proved in [7 Theorem 6.2]. 

Theorem 2.7 Let (Pe) £e (o g) be a family of normal states on J$f parametrized by £. Assume Tr[p £ N 5 ] < C§ 
uniformly w.r.t. £ £ (0,e) for some fixed 8 > andCg £ (0,+°°). Then for every sequence (£ n )n€N with 
lim„^oo £„ — there exists a subsequence (£ nk )keN cind a Borel probability measure jl on 3f such that 



limTr[p £ 



Wevh 



2? 



b{z) dn(z), 



forallb£U peP .^- 1 (Jt b (p2?)). 

Moreover this probability measure jj. satisfies / \z\ 2S djl(z) < °°. 

Definition 2.8 The set of Wigner measures associated with a family (Pe)ee(0,e) ( r ^P- a sequence (p £|I ) n ^) 
which satisfies the assumptions ofTheorem \2.7\ is denoted by 



3 £ ,ee (0,e)), [resp. J?{p En ,n eN)). 
Wigner measures are in practice identified via their characteristic functions according to the relation 

.,e £ (0,6)) = M & HmTr[p £ W(V2^)] = 

£^0 
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The expression ^(p £ ,e £ (0,e)) = {/J.} simply means that the family {ps)ee(0,e) i s "pure" in the sense 



limTr 

£^0 



for all cylindrical symbol b without extracting a subsequence. Actually the general case can be reduced to 
this after reducing the range of parameter to £ £ {£„ k ,k £ N} . 

A simple a priori estimate argument allows to extend the convergence to symbols which have a polyno- 
mial growth and to test to Wick quantized symbols with compact kernels belonging to (see Q 
Corollary 6. 14]). 

Proposition 2.9 Let (Pe) £e ^o £ ) be a family of normal states on jSf(Jf) parametrized by e such that 
Tr[p f N"] < C a holds uniformly with respect to £ £ (0, £) for all a G N and such that Ji(p £ ,£ £ (0,e)) = 
{jli}. Then the convergence 

]imTi\p e b^ mntized }= [ b(z)dp(z), (18) 

£^0 J& 

holds for the Weyl quantization of any b £ S v with p £ P and V £ [0, 1], and for the Wick quantization of 
any b £ &>™ lg (2r). 

Wigner measures are completely identified by testing with Weyl-quantized observable and possibly by re- 
stricting to some countable subset {J ne j$D„ tPn where D n>Pn is a countable dense subset of (y$b(p n 3?)), 
and (/?n)neN is a sequence of P such that sup„ GN p„ = Zy (see [7]). One may question whether testing on 
all the b wick with b £ ^ g (^) also identifies the Wi gner measures. When 3? is finite dimensional, this 
amounts to the well-known Hamburger moment problem of identifying a probability measure v on M from 
its moments a n = J M x"dv(x), n £N, for which uniqueness fails without growth conditions on the sequence 
(#«)neN (E3 0), which can be translated in our case to growth conditions of (sup eG ( ^ Tr [p £ N a ]) ae N- 
We shall circumvent this difficulty, by identifying the Wigner measures in two steps by approximating the 
states (pe)ee(0,s) by states {p £ pp ) £ £io, £ ) f° r which the growth condition is satisfied. We shall reconsider the 
moment problem later, but the comparison argument is given below. 

Proposition 2.10 Let (Pe)ee(0,e)> J ' = 1 , 2, be two families (or sequences) of normal states on Jff such that 
Tr p £ N s < Cg uniformly w.r.t. e £ (0,£) for some 8 > and Cs £ (0,+°°). Assume further <M (pi, £ £ 
(0,e)) = {nj} for j= 1,2. Then 

|Ml-M2|<limM|pi-pf|^i (jr) . 

£->0 v ' 

Proof. For a symbol b £ y cy i(3?), the finite dimensional Weyl semiclassical calculus says \b Weyl \^^,yf) < 

\\b\\^+ff b {e) with ||£||oo= \\b\\ L - {pjr) . This implies for a given b£5? C y l (%), 

| / b(z)d(^-ll 2 )(z)\ = lim Tr[(p] - p 2 £ )b w ^ 1 } < ||ft||-liminf |p] ~p 2 e \^ m ■ 

J %e £^0 £-5-0 v ' 

The measure p.\ — \ii is absolutely continuous with respect to the Borel probability measure ^' ^ 2 . Hence 
there exists a Borel function Aonf such that \i\ — \i 2 = X (z) ^"t^ 2 with the additional property | X (z) | < 2 
Ml t ^ 2 -almost everywhere. But for any Borel probability measure v on it was checked in [7| that 
<5*L/(iF) is dense in L p (3f,v) for p £ [1,°°) on the basis of a general measurable version of Stone- 
Weierstrass theorem (see for instance |fl9ll ). Hence there exists a sequence (j3 n )neN m -^cyl(^) sucn that 



X 

and after extraction lim^oo p\ k (z) = j^(z)i{x^o}(z), Ml ^^ 2 -almost everywhere. By setting b% = 2 
we get a sequence (b^keN sucn that 

Vk£N, b k £,Y cyl and ||A t ||-<l. 
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We conclude with 



J \p x -p % \= Jjk{z)\d 



Mi + P2 



(z) = 



Urn / ^(zJdOi 1 -/! 2 )® 



□ 

When the two sets ^{p ] £l e G (0, e)) have more than one element, the extraction of subsequences, (e n ) n€ ^, 
can be made simultaneously and the result has to be modified into: 



inf 

(m ,/i 2 )e-#(pi ,ee(o,e))x.4fr(p?,ee(o,e)) 



|Ml -M2I < Hmsuplpj -Pe\sfi(jtf) ■ 



(19) 



2.6 Wigner measures and the BBGKY hierarchy 

The compactness condition b G JS?) f° r the Wick quantization in Proposition 12.91 is not a technical 
restriction and the convergence is no more true for a general b G ^ a i s (3f). It was identified in [7| as a 
"dimensional defect of compactness" and illustrated with examples, one of them being related with the 
Bose-Einstein condensation of the free Bose gas. 

This terminology comes from the idea that this defect of compactness does not come from the infinity in 
the phase space like in the finite dimensional case (see [53 1[28|) but from the non compactness in the norm 
topology of balls in infinite dimension. Actually this was made more accurate in [8 1: under the assumptions 
^#(p £ ,e G (0,e)) = {p} and Tr [p E N k ] < l k , we proved (T) =>> (P) with 



(P) 
(T) 



V77 >0,3P T? gF 



limTr 

£^0 



Peb 



Wick 



b(z) dp(z); 



Tr[(l-r(P„))p e ]<T}, 



where (T) appears as a quantum Prokhorov condition (or tightness condition in the strong topology). 
The condition (P) which will be simplified in the next subsection, actually contains, for all a G N, the 
uniform bound w.r.t. e of Tr[p £ N a ] since N" = [(|z| 2 ) w ' ] a . It implies actually a strong relationship 
between the Wigner measure formulation and the convergence of reduced density matrices. 

Proposition 2.11 Assume that (Pe)ee(o,e) is a family of ^f l (Jif) with p £ > 0, Tr[p e ] = 1, 
(0,e)) = {p} with the condition (P) and assume p ^ 8q. Define for p G N 



,£ G 



7o 



Lp) 



y\z\ 2 Pdp(z) 



\z® p )(z® p \ dp(z) 



(20) 



Then for all p G N, the reduced density matrix Ye converges to yi p ' in the Jf^-norm. 

Proof. For p = 0, the result is nothing but \ = J p = lim £ ^oTr[p £ ] = 1. 
For p G N*, the condition (P) with p ^ 5q says first 



lim Tr 

£^0 



Pe{\z 



2p\Wick 



ftp 



2? 



dp(z) > 0. 



Hence, the reduced density matrix y £ is well defined according to Definition 12. 61 for £ < e p small enough 
(with Tr[p £ ] = 1). The condition (P) gives the general convergence: 



lim Tr 

£^0 



i p) h 



= lim 



Tr [p e b 



Wick' 



J* Hz) dp( z ) 

e " Tr \p e {\z\ 2 P) mck ] J r \z\ 2 " dp ( Z ) 



Tr 



for all b G ^p, p (i2°), where the last equality is a p -integration of the equality of continuous functions 



b(z) = (z® p , bz® p ) = Tr [\z® p ){z® p \b] 



This proves the weak convergence of Ye to Yq m ^ (V But since y( p ' E an d Yo are non negative 
with Tr[y £ p '] = 1 = Tr Yn , this implies the norm convergence according to |5O]|4]|2O|)0 



□ 



2 In a more general framework, it is said that J£ 1 (\J P 3*) has a uniform Kadec-Klee property (see [40] and references therein). 
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2.7 A simple criterion for the reliability of Wick observables 



The proof of Proposition ^. 1 1 l ean be adapted in order to make an equivalent condition to (P) with a weaker 
and easier to handle formulation: 



(PI) VaeN, limTr[p £ N a ] = / |z| 2a dfl( 

£^0 hi 



z) < +°°. 



Proposition 2.12 For a family (p e )ee(o,e) ' n ^£ l (^) such that p £ > 0, Tr[p £ ] = 1, ^#(p £ ,£ e (0,e)) = 
{jli}, the condition (P) and (PI) are equivalent: 



VaeN, limTr[p £ N a ] = / \z\ 2a dfi(z) ) Vfo e ^(Sf), HmTr p £ b Wick = b dp 

Proof. The condition (PI) is a particular case of (P). Let us prove (PI) =>■ (P) . 
We start with two remarks: 

• For/teN*, (|z| 2 *) Wfc * = N(N-£)...(N-(fc-l)e). Hence the condition (PI) is equivalent to 



Va e N, lim Tr 

£^0 



p e (\z\ 2a ) Wick 



z\ 2a dn(z) 



2? 



• For p — (resp g = 0) the operators in J?(C,\/ q 3C) (resp. in J?(\/ p JT,C)) are compact and 
@>0, q {&) = 9>%j\&) (resp. ^ P ,o(3f) = &~ fi (3?)). Hence the convergence lim £ ^ Tr[p £ fc w,c *] = 
fbdjj., is consequence of Proposition ^. 9| when p = or q = 0. 



According to Proposition ^. Ill there are two cases. 

If H = So: Then for b e ^ v ^(^), p e N*, such that > 0, the inequality < b < \b\&> Iyp % and the 
positivity ( fTTT i says 



< limTr 

£^0 



p £ £ 



< lim l^lp.pTr 



£^0 



PeCkl 2 ')™ = / \z\ 2p 5o(z)=0. 



For a general b e & PjP (3f), p e N*, the decomposition £> = £>k, + — t>R- + ibj t + — ibi- with all the b, > 
now gives 



Vp e N*,Vfr e 9 p A2r)i HmTr 

£— >0 



Wick 



= 0. 



For p q, p,q e N*, write 



Tr 



p £ /7 



Tr 



'pl l2 (pl' 2 b w H <Tr[p e p Tr 



p £ b Wick b Wick '* 



1/2 



Proposition El says that b Wick b Wick >* = lf i=Q ^b.dfb belongs to 0^^(^) with an <?(e) term in 
^oo(^)- We have proved 



Vp, ? eN*,V^e HmTr 

£^0 

while the cases (0,g) and (p,Q) are already known. 



Peb 



Wick 



= 0= / fe(z)5 (z), 



If jU ^ 5q: Then we know by Proposition ^. 1 ll that lim £ ^o II jf" 1 — 7q P ' ||^i = 0, which implies 



V/?e #„(f), limTr p £ £ 

£^0 



= lim Tr 

£^0 



f/'b =Xr y {p)~ b = J bQdnfc). 



Let us consider the general case b e 8P p ^(2?). The above convergence is still true when the kernel b is 
compact by Proposition ^. 91 Consider now a general b e S? PA (^)- Since / % \z® p ){z® q \ dji(z) is nuclear 
(or trace-class in V 9 & © V p ^0. f° r any n e N there exists a compact operator b„ e jf°°(\/ p V 9 
such that |fc„|^ M = \b n \j?(\jp 2?,\ji sr) = \b\jc(yp 2?,\ji 3t) = \ h \& PA and 



(b(z)-b n {z))dii(z) 



Tr 



\z® p ){z® q \dn(z){b-b n ] 



< 



n + 1 
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\r(p+q) 



2? 



The Lebesgue convergence theorem with 

V«GN, \b{z)-b n {z)\ r <(2\bW p J\z\ r{p+cl \ 

Vze2f, lim b n {z) = lim ( z ®« ,b n z® p ) = b(z) , 
n—^oo n — y°° 

yields 

lim / \b(z)-b n (z)\> -dn(z)=0. 
Set r\ r (n) = f^- \b(z) — b n (z)\' dji(z) and use again the Cauchy-Schwarz inequality 



Tr 



p £ {b wick -bf ck ) 



<Tr 



p £ {b Wlck -bT ck ){b 



1/2 



Owing to the result valid when p = g we deduce 



lim sup 



Tr 



p £ {b wick -bf ck ) 



< 



3f 



\b(z)~b n {z)\ 2 dn(z) 



1/2 



,1/2 



Since for n £ N fixed, lim £ ^ Tr [p £ bjf !t *] = J„ r b„(z) d/j.(z), we deduce 



Vn £ N, lim sup 

£^0 



Tr 



p £ £ 



b(z) dji(z) 



< 



n + l 



while the right-hand side goes to as n — > °° . 



□ 



2.8 States localized in a ball 

The condition, Tr [p £ N a ] < X a for all a £ N, used in (8) is actually equivalent to 

= l[o,A](N)p e l[ ,A](N) 

(locate the spectral measure of p £ for the self-adjoint operator N). Such an assumption remains an important 
step in the present analysis, and N = (|z| 2 ) w ' suggests that such a state is localized in ball of the phase- 
space. 

Definition 2.13 A family (p £ ) £6 (o j£ ) (or a sequence (p £n ) n en) of normal states on is said to be localized 
in the ball of radius R > 0, if p £ — l[o,«2](N)p £ l[ OS 2](N) for all e £ (0,e) . 

The meaning of the geometric intuition contained in the terminology "localized in a ball of radius can 
be made more accurate. 

Lemma 2.14 For a family '{Pe)ee(0,e) ( or a sequence (p £n ) n efi) of normal states on localized in a ball 
of radius R > 0, all its Wigner measures are supported in the ball {\z\ < R} . 

Proof. A family (Pe)ee(0,e) localized in a ball of radius R satisfies Tr [p £ N 5 ] < R 2S for all 8 > 0. Therefore 
the set of Wigner measures ^(p £ ,e £ (0, e)) is well defined and the convergence after extraction can be 
tested with Weyl-quantized cylindrical functions in the symbol class S v p introduced in ( fT~4T > for any p £ P. 
Let ji £ (p £ , £ £ (0,e)) be associated with the sequence (£„)„ g n- F° r any finite rank projection p £ P, the 
Wick quantized operator (Ipzj 2 )^^ is Np ®Ir s ((i-p)&) where N p is the number operator on T s (p2f) and 
equals (\z\ 2 ^ — C p e) Weyl in the finite dimensional framework of pfOf. For any cut-off function % £ ^°(R.) 
such that % = 1 on [0,^ 2 ] , the finite dimensional Weyl semiclassical calculus tells us (1 — #)(N P ) = (1 — 
X)(\z\ 2 pS -) WeyI + ^p(e) in Sf(r s (p3f)). Further the commutative decomposition N = fip'Si Ir s ((i-p)Sf) + 
h s { v 2f) ® > Np®/ rj ((i_p)jr) and choosing % decreasing on [0,+°°) implies 



(l-X)(\pz\ 2 ) Weyl + P (e) < (l-z)(Np®/r( (1 - p) .r)) < (1-*)(N). 



We deduce 



2\\Wevl 



< lim Tr 



0< / {\-x{\pz)Y)dll{z)= HmTr p £ (l - x(\pz\ 2 )) 

Hence the measure ji vanishes outside a cylinder {\pz\ > R}- This yields the result. 
With such localized states we can solve the moment problem. 



P £ „l[o,*2](N)(l-X(N)) 



= 0. 



□ 
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Proposition 2.15 Let (pe)ee(0,e) be a family (or a sequence (Pe„) n en) of normal states on localized in 
the ball of radius R > . If there exists a Borel measure jl on 3f such that 



Vfce^(iF) : 



limTr 

£^0 



Wick 



= / b{z)dn{z)., 

J Si- 



then 



^#(p £ ,ee(0,e)) = {^} 



Proof. Although this is shown in Q Proposition 6.15], we provide here a different proof. 
Let p £ P and consider the direct image by p of the measure pi : 

V£ g &{p2T) , \i P (E) = 1 r , (E) (z) dpi (z) , 

where &S(j>3T) denotes the Borel ff-set on p^ . 

For any b £ ^aigi-^)' sucn tnat b(pz) = b(z) we have 



lim Tr 

£^0 



Peb 



Wick 



pSf 



b(z) dp p (z) 



This holds in particular when b(z) = \pz\ 2k with b Wwk = N* + @{e) <N k + 0{e) with 



\z\ 2k dn P {z) < limTr 



< limTr 

e->0 



p £ N A 



<R 



2k 



Hence all the moments f \z\ 2k djX p (z) are bounded by R 2k and the finite dimensional moment problem 
applies (see Il47l ll5l): \i p is completely determined by the set of values | rfjUp , polynomialj . Let p! 

be a Wigner measure of the family (Pe)ee(o,g)- It is supported in the ball {z G 3f , \z\ < R} so that its direct 
image by p, fl' is supported in the ball {z £ p2f, \z\ < R}- Moreover there exists a sequence (e n ) ne ^, such 
that 



VbeS v p ^, 



lim Tr 



PeJ> 



Wcxl 



p T 



b(z) dv> p (z), 



where the b Weyl can be replaced by b Wlck for any polynomial b such that b(z) — b(pz) according to the finite 
dimensional comparison of the Weyl and Wick calculus in ( fTBI l. We deduce p p = \i' Since this holds for 
all the p £ P, this ends the proof. □ 

Let X be a continuous cut-off function supported in [0,1], with < % < 1 and such that % = 1 in [0, Aj. 
Within the assumptions of Theorem [27] and especially Tr[p £ N 5 ] < C§, the difference between the state p E 



and the localized state p £ 



— I — " — \%{iBi)p£%v5z) can be made arbitrarily small according to 



Vee(0,e), |p e -p, 



< 



(R/2)*-C S ' 



(21) 



where the right-hand side vanishes as 7? — ► «>. Then the comparison result in Proposition 12 . 1 01 or its variant 
( fl~9b says that the Wigner measures (Pe)ee(0,e) can be identified by its approximation by states localized in 
balls: 



inf 



Im-m'I< 



(R/2) 2S -C S 



(22) 



Then the question arises whether the family (p;r ) e6 (o,g), or an extracted subsequence, satisfies the condi- 
tion (PI) (or equivalently (P)) if the family (pe) ee (o,e) does. 

Proposition 2.16 Assume that the family (pe)ee(0,e) of normal states on Jfi? satisfies ^£(p £ ,e £ (0,e)) = 
{jj.} and the condition (PI) . Let the function f £ ^ ([0, +«>),R) be polynomially bounded such that the 
quantity Tr[p £ / 2 (N)] is uniformly bounded from below for £ £ (0, e) . Then the family (pi)ee(0,e) given by 

(0,8)) = { /fi($fdn } an d satisfies 



n f— l 

^ £ Tr[p E /2(N,j 

the condition (PI) 



/(N)p e /(N) has a unique Wigner measure 



13 



We will need the next lemma 



Lemma 2.17 Let the family (Pe)ee(Q,e) ( or a sequence (pe n )neH) of normal states be localized in the ball 
of radius R and assume the condition (PI) with ./#(p £ ,£ € (0,e)) = {/i}. Then the equality 



limTr 

£^0 



e a ^p e e a ^b wick 



J" 



b(z) dn(z) 



(23) 



holds for all a\ , tt 2 € C and all b € SP a l%(^) ■ 

Proof. The right-hand side of d23l is the sum of the double series 

^ (ai)*(c&)** 



ki\k 2 \ 



\2k l +2k 2 



b(z) dfi(z), 



for jA is a Borel probability measure supported in {|z| < R] and b is a polynomial function. 
Due to p e = p £ lr S 2](N), the sum 



JK 2 ,k - 

and the remainder term 

R K2 , k = Pee a ^(l+N) k 



k 2 =0 



k 2 \ 



SR 2 ,k — 



(i-0* 2 



K 2 \ 



p e (a 2 N) /f2+1 e a2,N (l+N)* dt 



satisfy 



1 [0,R2]{N)S K2M ^S K2 ^ k with \S K: 



and l [0R2] (N)R K2 . k = R Kl . k with \R Kl . k \^(.X>) < ^ i 1 + R 



3 2\k 



(\a 2 \R 2 ) K i +1 
(K 2 + l)\ 



Repeating the same estimate on the left hand side with Sk-, . k and R^ k instead of p £ implies that the Jzf 1 ( Jif) 
norm of 



(1+N)* 
is bounded by 

e \a 2 \R 2 +\a l \R 1 ^ + R 2-j2l< 



«,N n «,N & £ («lN)*' (« 2 N)^ 



kl =o k2 =o 



k 2 l 



(1+N)* 



(|ai|* 



2v£-,+l 



(l«2|/? 



l«l|^ 2 ) 



\a 2 \RT 2 



2\K ? + 1 



(#2 + 1) 



(^1 + 1)1(^2+1) 



kuWick 



which vanishes as min(Ki,K 2 ) — ► °° . We conclude with a 5/3-argument after noticing that (1+N) b 
(1 + N)~ A ' is bounded for k>kb and that the convergence ase->0 holds for € ^ a i g (3^) fixed and for 
the finite sums L^' =0 £jfe 2 =o owm g to the condition (PI) . □ 

Proof of Proposition [2JjD Let Cf > 1 be a constant such that Tr [p £ / 2 (N)] > ^ and sup iG ^ +00 )/(s)(l + 

s)~ v < C f . The inequalities Tr[p/N a ] < C^Tr [p £ N a (l +N) 2v ] , a G N, ensure that the family (p{)ee{0,e) 
admits Wigner measures without any way to identify them for the moment. So take a sequence (e„)„ G N> 
such that lim„^oo£„ = and ^#(p£,n 6 N) = {/lA . We first prove that the sequence (p{ n ) n eN satisfies 



/ 2 (kl 2 )M 



the condition (PI), then check that = — ; — - , 
then when / is compactly supported, and finally conclude with approximation arguments 



in the cases when (p E ) 



ee(0,e) 



is localized in a ball and 



1) The condition (PI) for the sequence: The uniform control of Tr 
/ 2" \z\ 2a djj.f(z) < +°° and the Proposition ^. 9| says that the convergence 



p£n« 



< C a , a 6 N, implies 



lim Tr 



Pib mck 



J" 



b(z) d\i s (z) 
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holds for any b G A3f) with a compact kernel. In particular for b(z) = \pz\ with /j£P and k € N, 



lim Tr 
while we assumed 



I (i\Pz\ 2 ) Wkk 



lim Tr 

H— S-oo 



/ (\\2k\Wick 



ST 



\pz\ 2k dn'iz), 



' P e„b Wkk 


-L 







b( Z ) dpi(z). 



HmTr 

n— >oo 

Fix a G N* and take 5 > 0. By Lebesgue's convergence, there exists p G P such that 

/ \\z\ 2a -\pz\ 2a \ diif(z)<8 

and / || z | 2a -|pz| 2o: |(l + |z| 2 ) 2v ^(z)<5. 

Jar 



(24) 



(25) 



Remember that (\pz\ ) = N P ® /r t ((i - P )iT) = N p with N" < N" where both sides commute with /(N) 
and we get: 



0<Tr 



p{(N a -N«) 



< C/Tr 



f(N) (N« - N«) 1 / 2 p £n (N a - N«) '/ 2 /(N) 



< 



C/|/(N)(1+N)- V k (jr) 



xTr 



(l+N) v (N a -N p ! ) 1 / 2 p £ „(N a -N ; 'J) 1 / 2 (l+N) v 



< C}Tr[p E „(N a -N£)(l+N) 2v ] . 
But we know by d25l > that the right-hand side converges as n — > °° to 

/^(|z| 2a - |/>z| 2a )(l + |z| 2 ) 2v dii( Z ) < C}8 , 

while (El with (|/Jz| 2 ) Wi ' c * = N p gives 



lim Tr 



\pz\ 2a dpf{z). 



Hence there exist njGN such that 



V« > rig . 



Tr 



■ir W 



2 « 



< (0 + 1)5. 



From Jj. ||z| 2a — |/?z| 2a | (z) < 5, we deduce 



lim sup I Tr 



p{n" 



5f 



z| za dn(z)\<(C} + 2)8. 



Letting 5 — » ends the proof of this part. 

2) Identification of p.? when (p e)ee(Q ,e) is localized in a ball: Assume that (pe) e e(0,§) i s localized in 
ball of radius R > 0. The Lemma l2.17l tells us 

Vfi , h e K, Vfc G &alg(&) , lim Tr [e" 2N p £ „e' YlN /7 W, ' £:A - 



while the uniform boundedness of (1 + N) *i>fc w,e *(l + N) A * entail 



Tr 



< C b Tr 



p £n (l+N) 2k » <C b (l+R z ) 



2\2k b 



Hence for / G (U (R)), we get 



lim 



Tr[/(N)p £ „/(N)^"*] _ /^/(|z| 2 ) 2 Kz)^M(z) 



Tr[/(N)p £ „/(N)] 



J^f(\z\ 2 ) 2 dn(z) • 
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We have proved: 



Vbe^ alg (2f), limTr p{i" 



f*f(\z\ 2 ) 2 b(z) dji{z) 
J T f(\z\ 2 ) 2 dn{z) 



The part 1) and l[o/;2](N)p£l[ 0/ j2](N) = p[ n ensure that (p£)„eN satisfies the sufficient conditions for 
solving the moment problem (Proposition ^. 151 ) and ^ 



:l 2 ) 2 / 



in this case. 



3) Identification of \i$ when / is compactly supported: Assume that / e ^ ([0,+°°)) is supported in 
[0,-Ro]. Consider for % e "^([O, +«>)), < % < 1, % = 1 on [0, 1/2] and for/? > 0, the truncated states 



Tr 



For > 2/?o, we have 



Vn 6 N* 



p{ = 



N. .N. 



Tr[p|,/ 2 (N 



r/(N)pf„/(N). 



By extracting a subsequence we can assume ^# (p§ n ^,k ENj = {^ R }, and Part 1) applied to (p R lk )keN 
ensures that the pair (p{ ,p| ) fulfills all the assumptions of Part 2) if / e ^([0,+°o)) n &~ l L l (M). 



/./ annolc /(N 



Thus the measure [i 1 equals j^^^^-j^ ■ From the comparison d22l we know / \[i R — /i \ = G(R 1 ) and 
since / is a bounded function 



l/(kl 2 )V 



c 

< — . 

- R 



J\f(\z\ 2 )\ 2 dLL 

Taking the limit as R — > gives the result when / € ([0, +<»)) PI J^ -1 /, 1 (IR) . Removing the condition 
f £ J^'L^R) is obtained by a comparison argument between p[ n and p{ l n with /f e ^? n ^" _1 L' (M) and 
sup se[0l+H 1/ W - /* (*) I < • for *e N ■ 

4) Final approximation argument and uniqueness of jx^: Consider now the complete problem with the 
extracted sequence (p£)„ G N- We again use the cut-off xijp) but now to compare p{ n with p{£ ■ After 

extracting a subsequence, we can assume ^# (Pe^ R \k ENj = i\l^ R 2 '| . The pair (p{^ R '\p£ n ) 
fulfills the assumptions of Part 3) and 

f 2 (\z\ 2 ) X 2 (R- 2 \z\ 2 )n 



fx(R- z -) 



J ' f 2 M 2 )x 2 {R- 2 \z\ 2 ) d» 



But from the inequalities f(s)(l ~x( R ~ 2s ))( l + s)~ v ~ l < CR 2 and Tr [p e (l +N) 2v+2 ] < C v we deduce 
the uniform estimate: 



Again the comparison argument d22l gives 



f _ fx( R - 

n k n k 



r' 

Zl 

^(30) ~ R 2 ' 



< 



! 



Jf 2 (\z\ 2 )x 2 (R- 2 \z\ 2 )dn 



< 



r' 
Zl 
R 2 ' 



and we take the limit as R 



tion (PI) was checked in Part 1). 



°°. We have proved ji s = jwrpw^ f° r an y sequence extracted from 

— $ ^ — — ^ while the condi- 

□ 



(p/) £G ( £) with a single Wigner measure. This proves ^(p £ , £ € (0, e)) = j ^[^j^ } 



3 Dynamical mean field limit 

Let Q be a real-valued polynomial in ^^(iF) given by 



<2=£e 7 -, with Qje&jji*)- 

7=2 
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We consider the many -body quantum Hamiltonian for a system of bosons 



H e = dT(A) + Q Wick , (26) 

with A a given self-adjoint operator on 3f. Here Q w,ck is the operator YIj=i Qj' ck with Q w,ck given by 
0. Clearly, H £ acts as a self-adjoint operator on the symmetric Fock space Jif. When 3f = L 2 (R d ), the 
Schrodinger Hamiltonian A = —A + V(x) and the semi-relativistic Hamiltonian A = \/—A + m 2 + V(x) are 
among the typical examples (e.g. ETTn . 

3.1 Existence of Wigner measure for all times. 

The first step to prove Theorem 1 1.1 1 is to show the existence of Wigner measures for all times. This is 



accomplished in the Proposition 13 . 3 1 by following the same lines as in the proof of Theorem l2.7l For this 
task two useful lemmas are stated below with the first one being proved in [7 Proposition 2. 10]. 

Lemma 3.1 For any b G 2? a i g {?£) we have: 

(i) b Wwk is a closable operator with the domain of its closure containing 

^ = vect{W(<p)y/,y/€^„,<p€ JT}. 

( ii) For any <p G 2? the identity 

W(S)*b wu *W(Z) = (b(z+^)) wick 

v2 

holds on Jf with b(- + ^fi) G &aig{2?) ■ 

Lemma 3.2 For any k G N there exists a e-independent constant C k > such that 

W^Y(N) k W^)<C k (E) k (^ k (N) k , (27) 
for any £, G 3f and uniformly in £ G (0, e). 

Proof. Since N is a self-adjoint operator, the functional calculus provides the inequality 

(N)*<(l+N)*. 

Therefore, it is enough to prove ( |27] | with (N) in the l.h.s replaced by (1 +N). The Wick calculus in 
Proposition|2~4ltell us that (1 + N)* is a Wick operator with symbol b k (z) in ® k - Q i.e.: 

h(z) = £ (z^'b^z^) with bf G 
j=o 

Now, applying Lemma lXTl vields 

w(sy(i+N) k w(S)=w(tyb¥"w(S) = {b k (z+^)) Wi *. 

A Taylor expansion of the symbol gives us 

bk (z+^) = t^ljD^b k (zm, 

v2 j=0 j\V2J 

with dW is the jth derivatives and D^b k (z)[£,] G ® k m i =a ^m,n{^)- So, by the number estimate ( TO we 
can derive the following bound 



(N)-*/2 (D^b k (zM]) W ' Ck (N)^/ 2 



<c k (^y 
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with Q only depending on k € N. Hence, we obtain 

1 (ie)J 



(N) 



-k/2 



i=0 jW2J v 



Wick 



(N) 



-k/2 



<C k {e) k {^)\ 



with Q only depending on k € N. Thus, we conclude that W(<^)* (1 +N)' : W(^) as a positive quadratic 



form is bounded by C k (e) k (N) k 



□ 



Proposition 3.3 Let (Pe) e g(o g) ^ e a family of normal states on Jtf' satisfying the uniform estimate Tr[p e N r ] 
< C r for some r > 0. 

Then for any sequence (£„)„ g n !n (0> *0 f/zaf lim„_i.oo £„ = f/iere exists a subsequence (£ nk )keN an d a 
family ofBorel probability measures {Ht)teWL satisfying 

e-^ H *' p £i /i H ^ ,n € N) = {Mr} , 



for any t € R. Moreover, we have 



Proof. We set 



/ \z\ 2r d^{z)<C r . 



p £ (t) = Pe e iL e H ° and p e {t) = e 'W) g -.-£ff £ p jLH e e -i|dr(A) 

(i) Consider for e > the function 

G e (f,D=Tr[p e (f)W(V2*£) 

Write for any (5, | ) , (f , 77 ) e R x 3? 
\G e (t,T])-G e (s,$)\ < \Tr[(p e (t)-p e (s))W(V27CT}) 
By differentiation, we get 



Tr 



peO) (w(\/27r77)-W(\/27T^) 



Tr 



(p e (/)-p e (s)]W(V2^Tj) 



< 



Tr 



p e (r')[G^,W(V2^T/)] 



with 2r' (z) = Q{e "' A z), while the second term is estimated by 



p e (s)(w(V27t7])-W(V27t%j) <(l+Cr) [W(V27trj)-W(V27t^)](N+l)- 1 



Tr 



Now, we claim that there exists a constant c > such that the r.h.s of (1281 1 is bounded by 



^-*l(Illfi;ll)Ie i - 1 |T?|'. 

j=2 (=1 

This can be proved by first writing 

Tr[p e (/)[e^,W(^Tj)]" - 

Tr [(N) r p £ (f')(N)'' ((N)- r W(V2^Tj)(N>'') {N)- r [W [V2%i])* Qf ck W {V2%r)) - Qf ck ] (N)" 
and second estimating the r.h.s of (|3TT i using Lemma [372] and Lemma IXTR ii) so that 



(28) 



(29) 



(30) 



(31) 



Tr 



p e {t')[Qy k \W{^/2m])] 



dt' 



<c\t — s\ sup 



l£ 



(N)- r [Q t ,(. + ^r,) W,ck ~Q7 ck m- r 
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Thus, the bound ( f30b follows from the number estimate in Proposition [Z 
We recall the inequality proved in [7 Lemma 3.1], 

[W(\/27rT7)-W(\/27z;£)](N + i;r 1/2 <C |tj-§| [min(e|Tj|,e|§|) + max(l,e) 

This leads to the following bound on the r.h.s of ( 



c<e>|T7-§| (i + Vl^l 2 

Thus, we conclude that V(i,|),(f, 77) eMxf, 



|G £ ( f ,T7)-G £ (,^)|<c^-,|(|77| + l) 2r +|77-^l V / |T?| 2 + |^| 2 J, (32) 

uniformly w.r.t. e € (0,e). Recall also that we have the uniform estimate |G e (s,<!;)| < 1. 
Now, we apply an Ascoli type argument: 

• Since M x is separable, it admits a countable dense set jV = {(t(, <jj^), £ 6 N}. For any I £ N the 
set {G e (t(, <^)} £G (o,£) remains in {a £ C, | O" | < 1}. Hence for any sequence (£„)„ g n such that e„^0 
there exists by a diagonal extraction process a subsequence, still denoted by (e„)„ e N, such that for all 
I e N, G £ „(^,^) converges in {a £ C,|a| < 1} as n °°. Set 

Gfe,&):= lim G e „fo,&) 

for all f el 

• The uniform estimate ( f32b implies that the limit G is uniformly continuous on any set 

^n{(t,z) eix jT: |f| + |z| </?}. 

Hence it admits a continuous extension still denoted G in (M x iF, | | Kx oj). An "epsilon/3"-argument 
shows that for any (t, G R x iF, lim„^oo G £(1 (f , ^ ) exists and equals G(t,£,). 

Finally for any t € R, G(t, .) is a norm continuous normalized function of positive type since 
G(t,0) = limTr[p e (?)] = l 

The positivity in the last statement follows by Weyl commutation relations ((6). Therefore, according to the 
Bochner theorem (e.g. |[T2l Corollary 1.4.2]) for any t £ R, G(f,.) is a characteristic function of a weak 
distribution or equivalently a cylindrical measure p., on 5° (see ll5D and also Q Section 6] for specific 
information). 

(ii) The fact that fi t are Borel probability measures satisfying 

£r(k| 2r )<C,.<~ (33) 



follows directly by 11461 Theorem 2.5 Chap. VI] or by part (iv) in the proof of Q Theorem 6.2]. 

(iii) Using ( TT3l we see that for any b £ based on a finite dimensional subspace p2? with p £ P 

lim Tr p En {t)b Weyl = Km / G^) J?[b](i;) L p (d!;) 

G(t,S)&\b](S)hVS)= I Hz)dfL t (z). 

P ar Jar 

Therefore, according to Definition 12. 8 1 of Wigner measures we conclude that 

VreK, ^f(p £ „(f),«eN) = {£,}. 

(iv) Finally the family of measures ji, which satisfy the claimed statement in the proposition are the push- 
forward measures 

Ht = (e~ ,M )*£r- 

Furthermore, an analogue of (l33l can be easily checked for the measures ji t . □ 
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3.2 Polynomial approximations of the classical flow. 

With the classical hamiltonian 

h{z) = {z,Az)+Q(z) = (z,Az) + £ (z® J , Qjz® J ), z € 9(A) . 

;=2 



the related nonlinear field equation is 



id t Zt =Azt + d z Q(z t ) 

Zt=0 = ZO ■ 



Actually this Cauchy problem is better studied when reformulated as an integral equation 

z, = e-' aA z - i f e- i( t-^ A d ? Q(z s ) ds, for z € S , (34) 
Jo 

which admits a classical ^-flow F, : K x — > Sf; 1) since the Qj are bounded a fixed point argument 
gives the local in time existence and uniqueness; 2) then the conservation \zt \ = \zq\ ensures the global in 
time result. As a classical ^°-flow, F is a ^-map satisfying F, +V (z) = F, o F s (z) and F,(z) solves ( f34b for 
any zef . 

Moreover, if Zt solves p4|) . and Q t (z) = Q{e~" A z), then w t — e" A zi solves the differential equation 

d ^ , . 

— w, = -«^e f (w f ). 

Therefore for any b € &pq(3f), the following identity holds 

-rb{w t ) = d l b(w t )[-id- z Q t {w t )] + d z b{w t ){-id- z Q t {w t )] 
dt 

= i{Q t ,b}(w t ). 
Hence, we obtain the Duhamel formula 

b(z t )= b,(z) + i jf {&, A} {e itlA z h ) dh . (35) 
A simple iteration in (l35l l. using 

{ft, , MK ) = {Qn , bt}(w ) + i P {Q h ,{Q, V Ml K )dt 2 , 

Jo 

yields 

b( Zt ) = b t (z) + i f {Q n ,bt}(z) dt! + i 2 fdh f dt 2 {Q h , 
Jo Jo Jo 

Therefore, by induction and after setting F f (z) = z ( , we obtain for any K > 1: 

6oF,(z) - Mz)+E*'* f dh--- P~ l dt k {Q tv {...,{Q h ,b t ] ■■■}}& 

j~! JO JO 



c fdh- f K ~ l dt K {Q tK ,{...,{Q tl ,b t }-}}(e itKA z tK )- 
Jo Jo 

With the polynomial Q we associate the norm 

II6II = max \Qj\# = max 12/1^^^ (36) 

and we note that \\Q t \\ = \\Q\\ for all t E M . Notice that the flow F, preserves the norm 

VzeiT, |F,(z)| = |z|, 
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and is gauge invariant 

Vz G X , V0 G R, F, (e i0 z) = e'% (z) . 

But for a given polynomial b(z), the map z i-> £>(zr) does not remain a polynomial. Starting from a polyno- 
mial b(z) G & p ,q(3f), we study polynomial approximations of b(zt)- 
Consider the expression 

b K (t,z) = b t (z)+£i k fdh-- f k ~ l dt k {Q tk ,{... 7 {Q h ,b t }...}}(z)= K £b k (t,z) (37) 

k=l Jo J0 k=0 

R K {t,z) = f fdh--- F~ l dt K {Q tK , b t }...}}(e^ A z, K ). (38) 

Jo Jo 

The two approximation results that we will use are given in the two next propositions. 
Proposition 3.4 For b G the polynomial b K (t,z) = Lf=o' bk(t,z) defined in ( f37b belongs to 

^j+pj+qi^) with the estimates 

\h(t,z)\<2^(p + q)(4P) k ||e||* \byjt\ k (z) 2k{r - 1)+p+9 - (39) 
Moreover, we have for R K (t,z) the estimates 

\R K (t,z)\<2^(p + q)(4rY \\Q\\ K \byjt\ K (z) 2K ^ + " + " . (40) 
Proof. With b G 9>p A {2?) and Q, = Y! j=2 Qj.t, the polynomial 

b k {t) = (if J* dti- dt k {Q tv {. . . , {Q n ,b,}- -}}(z) 
is the sum of (r— 1)* < r k monomials 

**(*) = L b kA f ) 

ae{2,...,r} k 

with b Ka (t) = (l) k J dh- J"' ' dt k {Q ak fk , { . . . , {Q a{ A , b t } . . .}} e&\ a\-k+p,\a\-k+q(^) ■ 
A consequence of Proposition ^. 4l savs for c G HPpirf(3f) , 

\{Q ai , h ,c}(z)\<r(p> + q >) \Q ai \^ \cU,, ( z f+*+m-i). 



(zV 

We deduce 



Mw)l < / o '*.-jf""'<*.' i (p+9)-(p+9+2t('-i))IIQII'W*„{z) M « +2| '" | - a 

< o + ,m- or | ( |. F £fe±i+2_- T i TTller IH ,„ 

with a = 2{^-\) anc ^ ^ denotes the Gamma function. Now, we notice the relation with the Beta function 

which yields d39l . 
The remainder 

R K (t,z) = i{Qr,b K } = F £dti ■ • • j'^ dt K {Q, K ,{- ■ ■ > {QnM ■ ■ ■}}(e i ' KA Z tK ) 

is analyzed like the term b k (t) . □ 

Proposition 3.5 Let pi be a positive Borel measure on 3? supported in the ball {|z| < R}, R > 0, then for 
any polynomial b G £P p ,q(3?), 

jR K (t,z)\ dpL{z) < (R) p+ «2W^{p + q )\by M [4r 3 ||e|| <fi^\t\]* . 
Proof. It easily follows from (140b . □ 
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3.3 Transport for a state localized in a ball 

The previous approximation result allows to prove partly Theorem 1 1.1 1 for states localized in a ball, intro- 
duced according to Definition 12 . 1 3 1 and studied in Subsection l2.8l 



Proposition 3.6 Let (p £ „),2eN be a sequence of normal states on ^ localized in a ball with radius R > 
and such that 



Vr G \—T,T] , J?{e~ l * H °»pe n e l * Ht » ,n G N) = {/*,} , 
and VcteN, lim Tr[p e „ N a ] = ( \z\ 2a dHo(z) ■ 



Then for all t € [-T, T], the probability measure jl t is the push-forward by the flow F r of the measure Hq, 
i.e., [l t = {f?t)*lki- Moreover the identity 



lim Tr 



~ i -k He «p en e i -k H£ "b quanti 



ZX'il 



b{z)d i i t {z)= / b{¥ t {z))d^(z), 

2 Js? 



holds for Weyl quantized cylindrical functions b G UpeP & {-^b (P^?)) an d general Wick quantized poly- 
nomials b G ^ a l g {3f). 



Proof. We set 



p £ „(f) :— e E « * 'e E " "p £ <? E « e "e e » v '. 



It is worth noticing that for all ret, the sequence (p £(1 (f ))„eN is localized in the ball with radius R. 
For a fixed G ^ pq (3f), differentiating with respect to t the quantity Tr[p £ (f )fr ], we obtain 



Tr[p £(l ( f )^]=Tr[p £ „(0)fr 



and replacing fr by fr r we end up with 



Wicki 



£n JO 



- Tr 



Tr[p £ „ (0 b"«*\ = Tr[p £ „ (0) b»«*\ + i / Tr p £ „ (,) {ft, fr f } 

Jo 



U7, /,' 



r e^ 1 

-Z 7 



Tr 



wCO ({&A} ( 



11 7, 7, 



(41) 



(42) 



Consider now the case when fr G &p q (Sf) with a compact kernel, fr G J£°°(y p 3f;\] q 2?) . Then we 
know that the left-hand side converges to fj»fr(z) djx t {z). The number estimate of Proposition 12.31 with 
Tr [N a p £ J < /? 2a implies that the last term of the right-hand side converges to as n — » °°. Finally the first 
term of the right-hand side converges to f ^ b(z) dflo(z), even when fr is not compact. 
We conclude that the limit of the second term of the r.h.s exists with 



b(z)dMz)= b t (z)dHo(z) + lim i Tr p En {s){Q s ,b t } 

3f Jsr n ^°° Jo 



Wick 



ds , 



and this initiates our induction process. 

Given K > 1, take the approximation b K (t) = Y,f=o bk(j) to fr(F r (z)) given in d37l i. and assume 



b(z)dii t (z) = / b K (t,z)dno(z) 

3f Jar 



lim/ fi / dty- / dt K Tr 
Jo Jo 



fr-i 



fe(fe)({G^ ! ---{<2 tl ,^}---}) mc * 

A simple differentiation with respect to tg_ gives for © € ZP a i g {3f), 



Tr 



U7,A 



Tr 



p £ „(O)0^ 



/ / flffc+iTr 
o 



' E — / Tr 



Pen(tK+l)({Qt K+1 ,®}) 
Wick 



Wick 



({Q, K+l ,®} 



(./) 



dtK-t 



(43) 
(44) 
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Hence, choosing © = {Qt K ,- ■■{Qt l ,b t }---} yields 



3T 



b{z) dnt{z) 



sr 



b K (t,z) dno{z) 



lim < i 



:K+\ 



•X 



du--- 



tK-1 



dtK Tr 



Wick 



;K+l 



dU ■ ■ 



'K 



dh--- I dt K+ {Yx 



Pe n (tK+i)({Q tK+1 ,---{Q tl M---}) 



Wick 



dt K+i Tr 



Pe n (tK+i)({Q tK+1 ,---{Q tl M---}) 



Wick 



1+ lim (II + III + IV). 



For any K, when n — > °°, the second term (II) converges to /«y®(z) djio(z) because the initial states p £ll (0) = 
p £n satisfies lim £ ^oTr [penC 1 *"*] = f^ c(z) dj^(z) according to Proposition ^. 121 Moreover, the third term 
(III) vanishes, when n oo, thanks to the number estimate in Proposition ^ . 3 1 and the fact that Tr [p£„N a ] < 
R 2a . Therefore, we have 



/ b(z)dfk(z) 



(t,z) dpL Q (z) 



lim i 



■K+l 



<K 



dt\--- / dtK+i Tr 



Pe n (tK+i)({Q tK+l ,--{Q h ,b t }--}) 



Wick 



By Proposition l3 .5 l and the fact that /io is supported in {\z\ < R}, we deduce 



3T 



b(z)dnt{z)- / b(F,(z))dn 



2? 



< (R)P^2W^(p + q)\by i , a \4r 3 \\Q\\{R)*-V\t\ 

fs~> ; 1 isWick 



(45) 





r'K-i 




lim / dt i • • 


■ / dtK Tr 


Pe n (tK)({Qt K 




Jo 



The number estimate of Proposition ^. 31 with the inequality ( f39T > of Proposition l3.4l implies 



(N) 



2 



({G fe ,--{G ?1 ,M--})^ (N) 



p+g(r-l) 



< 



2^( P + ? )(4rY||er 



This provides for the last term in the r.h.s of d45b the upper bound 



R) 



2*^b + ?)(4r 3 f ||Gf |*| 



(^lieiK*) 



7— [- with 5 < 1, taking the limit as K —> °° now gives 



For small times, |f| < T$ 



But according to Proposition 13.61 the measure \i t is a Borel probability measure supported in the ball 
{\z\ < R} which is weakly compact. Meanwhile cylindrical polynomials which are contained in 
because they are associated with finite rank kernels, make a dense set in the ( if (Z?(0,7?) vl , eflJ t,C) and there- 
fore in L (3f,djJ.) . Thus, we have proved 

Vte[-r s> r a ], p, = (F r )*po. 

Finally, since |F/(z)| = \z\ and [H e ,N] = 0, the pair ((p £(l (f))„ e pj,p f ) satisfies the same assumptions as 
((PgJneNjMo)- Since the time Tg depends only on Q and R the result extends to all f € M. □ 



3.4 Proof of the main result 

Gathering all the information of Section|2]and|3] we are now in position to prove Theorem ll.il 
Proof of Theorem ll.il 
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Let (pe)ee(o.e) be a family of normal states satisfying hypothesis of Theorem l 1.11 and let % € ^ ([0,oo),IR) 
be a continuous cutoff function such that < % < 1> X ( x ) — 1 iix< 1/2 and % (x) — if x > I. For R > 0, 
consider the family of normal states 



Pe = 



Tr{ X (N/R2)p £ x(N/Ri)} ' 
localized in the ball of radius R. By Proposition ^. 161 we know that 

,2[\„\2 /o2 



(h) VaeN, limTr[p £ *N a ]= / |z| 2a ^oW- 



Next, we use the notations 

Pe(t) ■■ 



p £ e 



and p*(f) =e~ i i He p§ e^ He 



For any sequence (£„) n eN there exists by Proposition 13.31 a subsequence (£„ k )keN an d a family of Borel 
probability measures (juf)f G K such that 

(0' ^(p4(0,*eN) = {M/ f } 

Mai 



(if)' VaeN, limTr[p* N a ] = 



i2a j,,fi 



2" 



d^(z). 



Applying now Proposition l3.6l with (/)' — (z'z)', we obtain that 

jets* (t),kzm) = {(F t )*n«}, 



(46) 



for any time t £ K. Since for any sequence (£„)„ g n we can extract a subsequence (Cn^eN sucn that (l46l > 
holds we conclude that 



(0,ce(0,e)) = {(F f ),Mo} i 



(47) 



for any R > and f e R. Again applying Proposition 13.31 for (p£) £ e(o,£)> there exists for any sequence 
(£n)neN a subsequence (£ nk )keN and a family of Borel probability measures (pi)teR such that 

Jt{ Penk {t),keH) = { i i t }. 

The identification of the measures {Ht)teVL follows by a 5/3 argument. For any b € ^ cy i(3f) based in pS 
p e P, we write 



Tr[p £ „ (t)b We *}- / ft(z)d(F,).f«, 



r 



< 



Tr[p % (0^]-Tr[pl(0^] 



Tr[p* (0^]- / b(z)dn' t 



2f 



[ b{F,{z))d^-j b(F,(z)W c 



(48) 
(49) 
(50) 



Each term (|48T> - (f50b can be made arbitrarily small by choosing R and k large enough and respectively 
using the bound d2TI ). the relation ( f4Tb and the dominated convergence theorem. So, we conclude that 
lit = (F,)*Hq and hence we have proved 

^(p £ ,ee(0,e)) = {(F r )*Mo}. 
Finally, the use of Proposition ^. 12"l with p £ (f) yields 



limTr 

£^0 



Pe(t) 



■Wick 



b oF t(z) dno(z), 



since lim £ ^ Tr[p £ (?)N a ] = lim £ ^ Tr[p £ N a ] = \z\ 2cc dHo = $ z \z\ 2ct dpL u for all a € N. The reformula- 
tion of this result in terms of BBGKY hierarchy of reduced matrices is a consequence of Proposition ^. 1 11 

□ 
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3.5 Additional results 

Although it was not written in Theorem 11.11 remember that the existence of Wigner measures contains a 
result for Weyl observables. 



Corollary 3.7 Let (p £ ) £e (o,e) be a family of normal states on 3ff satisfying the hypothesis ofTheorem \l.l\ 
The limit 

lim Trfe - '** p £ e l e H * b Weyl ] = f boFfz) dpo 

Jar 

holds for any b in the cylindrical Schwartz space ^^{Sf) , any t G M. and any b G S v V G [0, 1], p G P. 

The next result, shows that the class of observables can be extended to functions of Wick-quantized sym- 
bols. 



Corollary 3.8 Let (p £ ) £e (o,e) be a family of normal states on ^ satisfying the hypothesis of ' Theorem U .71 
Then 

i) The limit 

limTr[ e -'^ p £ e iL e H ° f(b wkk )] = ( f(boF t (z)) d^ (51) 

£^0 JSC 

holds for any f G (^t>(M.)) and any b G ^ Pj p(i2°) such that b* = b. 

ii) If additionally (p £ ) £ e(o,e) ' s a family of localized states on a ball of radius R > 0, then the limit ( 1511 ) 

holds for any entire function f(x) — YZ=0 a k xk over C and any b G 3 g p :P (3f) such that b* = b. 

Proof, i) Let % G ^([OjOo^R) be a continuous cutoff function such that < % < 1, %(x) = 1 if x < 1/2 
andj(x) = if jc > 1. Consider the family (p £ (t) = e~'^ He p £ e'i He ) £G (o, £ ) with 

„ z(N/^ 2 )p £ z(N/j? 2 ) 
P£ Tr[ Z (N/tf 2 )p £Z (N/fl 2 )] ' 

Let b G ^^(5°) such that £>* = fr, then fr 1 ^ extends to a self-adjoint operator on Jtf satisfying [N, fr^ 1 '*] = 
0. We claim that 

V6 G R, Tr[p£(t)e mWick ] = £ ^- 0*Tr[p*(f ) [b wick f] . (52) 

Thanks to the estimate 

Tr[p £ fl (f)(Z? w,c 'Y] = Tr[(N) M/2 p £ fi (r)(N) p * :/2 ((N)- p/2 fe w '''*(N)- p/2 ) ,; ] 

the l.h.s of d52i l is an absolutely convergent series uniformly in e G (0,e). Moreover, on can easily show 
the strong limit 

s -l im ttl 1 M (N) = e mW ' Ck 1 M (N) • 

Therefore, we see that 

£ ^e k Tr\pf(t)(b^ k ) k ] = f ^ fe Tr[pf(0(^)*l [0 ^(N)] = Tr[pf (r),^^] . 
This proves (l52l l and again by the uniform estimate d53l ) with respect to e G (0, e), we obtain 
limTr[pf(0 e ' e n = £^V KF f ( Z ))^Mo= fe^^d^. 
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Now, a similar 8 /3 argument as in the proof of Theorem ll.il 



< 



2 



jeb(F t (z)) 



>3f J3f 

using the bound 1211 . the relation ( |47] i and the dominated convergence theorem, yields the limit 



limTr[p £ (t)e 

£^Q k 



idb n 



jeb(F,(z)) 



2T 



By integrating with respect to ,^(f) € ^,(R), we end the proof, 
ii) The proof is similar to (i). Indeed, one shows 

TT\p e (t)f(b mck )} = f i a k TT\p e (t)(b 

k=0 



Wick\k] 



(54) 



with a l.h.s absolutely convergent series uniformly in e 6 (0, e). Letting e — » in (|54| ) yields the result. □ 

4 Examples 

We review a series of examples. Firstly, the propagation of coherent states and Hermite states is recalled. 
Secondly, bounded interactions occur naturally within the modelling of rapidly rotating Bose-Einstein con- 
densates, owing to some hypercontractivity property. Thirdly, the tensor decomposition of the Fock space 
allows to specify some Wigner measures for which the propagation cannot be translated in terms of the 
reduced density matrices without writing all the BBGKY hierarchy. Finally, the result of Theorem 11.11 
provides a new way to consider the Hartree-von Neumann limit in the mean field regime. 



4.1 Coherent and Hermite states 



The coherent states on the Fock space, T s (2f) are given by E(<^ ) =W '(■^^)f2 = e 



£2, where £2 is the 



vacuum vector of T s (2f), % £ 2? and [a(f),a*(g)\ =e(f,g)I. The Hepp method ([36 1[30| [31 1) consists 
in studying the propagation of squeezed coherent states a slightly larger class which includes covariance 
deformations. The normal state made with £■(<!;) is 

Pe G):=\wC§S)n)(w£§Z)n\. 

We proved in [7| that ^#(p £ (^),£ <E (0, e)) = {S^} and a simple computation shows that the property (PI) 
is satisfied: 

limJr\p e N k ] = \l;\ 2k = Ss(\z\ 2k ). 

A second example is given by Hermite states, also well studied within the propagation of chaos technique 
or other works (e.g., [44][ 13|[23|). They are given by 



(55) 



with <p £ 3f, \tp\%- = 1 and discrete values for £ = jj- We know from [7| that ^(p N ((p),N £ N) = 
{ Ik So" *V 8 <p ^9} where the rotation invariance is the phase-space translation of the gauge invariance of 
the Hermite states <p — > e' e (p . One easily checks the property (PI): 



lira Tr[ PN (<p)N k ] = 1 = — / \z\ Zk 8 e , e ^z) de 



' 2n Jo 

It is convenient to introduce a notation for this Wigner measure. 



26 



Definition 4.1 For <p € 3f, f/ie symbol 8^ denotes the Borel probability measure 

Theorem l 1.11 applies and the Wigner measures associated with 

(e-^pe^y^)^^ and {e- iL ° H *p N {<p)e iL e H *) e=mNeW 

are respectively 5^ and 5£ , where ^ f or <p, evolves according to the classical flow. 
For example when 

H e =dT(— A) + - / V(x — y)a* (x)a* (y)a(x)a(y) dxdy 

2 Jg2d 

with i2° = L 2 (M d ) the classical flow is the Hartree equation 

id t y= -Ay+(V* |v / | 2 )V / - 

We conclude by noticing that for such states (Pn(<p) and p £ (| )) the asymptotic one particle reduced density 
matrix )g (f) solves the equation 



-A + (V*nm),^ 
'o 



(56) 



y^(, = 0) = |<^l forp e (0, (resp. f l >(t = 0) = \cp)(cp\ for PAr (p) 
withn (1) (jf) = y^ ) 

'0 

4.2 LLL-mean field dynamics for rapidly rotating Bose-Einstein condensates 

The case of bounded interaction terms occurs exactly in the modelling of rapidly rotating Bose-Einstein 
condensates in the Lowest-Landau-Level (LLL) regime. The (LLL) one particle states can be described 
(see [2]) within the Bargmann space 



where L(dC,\) is the Lebesgue measure on C, h > is a parameter which is small in the rapid rotation 
regime and where the norm on 3f is given by 

I A = / C l/ Kl ,iv¥^g|i = -L|„ gl , „({,, ■ 

The multiparticle bosonic problem has been considered in BP and the (LLL)-model has been justified 
for the stationary states of such a system not only in the mean field asymptotics. The ^-particle states are 
elements of 

\J&=[f €L 2 (Cl,e-^ L(dQ)), d ( F = 0, ...,£„(*)) =F, Vcr e & k 



with the norm 



,2 L(dQ 



\C 2 



With or without the symmetry condition, ® 2£ and V * 2£ are closed subspaces of L (C|,e ~i<~ L(dQ) 
and they are the image of the orthogonal projection (add the symmetry for V 3f) 
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Within the modelling of rapidly rotating Bose-Einstein condensates, the one particle kinetic energy term is 
A = h^\d^ x and it is associated with 

0<E ki „(f) = (f,hC l d Q f), T . 
The standard one particle nonlinear energy is given by 



a / \u\ 4 L(dCi), a(Ci)=/(CiK 



c 

where a > is another parameter provided by the physics (see||2]), but more general energies can be 
considered ^ 

E NL (f) = t a p[ H 2p iKi), «(Ci)=/(CiK^K CC p >0. (57) 

p=2 JC 

The mean field Hamiltonian is thus given by 

h(f)=E kin (f)+E NL (f) = </,/!&%/> + t a p [ |/(Ci)| 2 '^^ L(rfCi). 

An important property of these nonlinear energies comes from the hypercontractivity of the semigroup 
£, ^ Q p rove( j i n ||T6l which can be written as 

\U\ L 2 P <C pXd \U\ L2 if U(Q=F(Qe-^, Fe®2f, pe[2,-H»]. (58) 

This implies that the nonlinear energy is a norm continuous polynomial with respect to / G 5° and therefore 
the nonlinear mean field equation 

idtf = E POpnJdal^^Hjni/ (59) 

p=2 



defines a nonlinear flow on the phase-space 3? according to Subection l3.2l (we refer the reader to [45 1 for a 
more detailed analysis of the nonlinear dynamics of the LLL-model) . 

Let us consider the second quantized version H £ of the energy h in T s (3f). The kinetic energy is nothing 
butdr(A): 

k 

dT(A) | v * x = e £ fc£ y d ? . = e/zC-^ • 



y'=i 



and the quantum Hamiltonian // £ is then 



H £ =AY{A)+Y j a p Qf ck (60) 

with Q P (f)=f i«(Ci)i 2p iK0= / i/wiv^iKonr,^). (6D 

Jc Jc 

The operator g /7 is easily identified after removing the center of mass in multiple integrals (see ATI for 
details) as 



Q P F(Q=n h p 



.7=1 



with £j = £j — + p + ^ p . One easily checks as well, by using additionally the hypercontractivity estimate 
dSD with p = +oo, that Q p e ^(\J P St). 

The propagation result of Theorem 11.11 applies for such a model for all initial states which fulfill its 
assumptions (boundedness of all moments and condition (PI)). 



28 



4.3 Fock tensorization 



We have already used, and it is the basis of the introduction of cylindrical observables, the fact that 

x 

T s {3f) ~ T S (3F\) ® r s (^) when = iFi © . The definition of Wigner measures introduced via cylin- 
drical observables, yields the next result. 

Lemma 4.2 Assume 2? = 3f\®3f2 an d let (pJ) ee (o,g)> (P^eeiOe) ^ e two families of normal states on 
T S (2Z\) and T s {^2) such that Tr [pjNf] < C§ holds uniformly for some 8 > and ^#(p £ , £ G (0,e)) = 
{/^} for £ = 1,2. Lef p £ be the state on T s (3f) identified with p\ ® p| in the decomposition T s (3?) ~ 
T s (3?i)®T s (3?2). Then the family (pe)ee(0,e) admits the unique Wigner measure jj. = jX l x jii 2 on the phase 
space 2? = 3?\ x . 

Before giving applications and variations on this result it is worth to notice that the identification of 
the "tensor" state p £ requires some care. It is not equal in general to p\ ® p £ since such a states does not 
preserve the symmetric Fock space T s (2?) . 

Here is a simple example, take (pi G and ^£^2 with \<pt\s- t = l,N\ € N, and set p e = \<pf e )(<pf 1 \ 
for i = 1 , 2 . The tensor states p 1 ® p 2 is the pure state | (pf Nl ® <pf N2 ) ( <pf Ni <g> ipf* 2 | in T, ( ) ® T s ( . 
It suffices to identify the vector q) v ( N i' N 2) g r i (^°) associated with . It is the symmetric vector 

in \/ Nl+N2 £f made with A^i -times <f>i andA^-times (fe and we can summarize the situation with 

1 Nf times 

<pf Nt = a*(<p e )...a*(<p e ) \a t ) mr s ($? e ), *=1,2, 

m v(NiJh) / (jVi+jV 2 )! f m ®*iffl m »% 

* = Y e(Ari+AW *i w 2 ! 

I iVl times times 

= ^=^===«*(<Pi)---«*(<PiK(<P2)---«*(<P2)|n) inr. 5 (^). 

The tensor decomposition is especially useful when iF is endowed with a Hilbert basis (e ,•) ; e j$* . An Hilbert 
basis of r s (3T) is (e va )a GU 7 =0 (N*y given by: 



with a natural multi-index notation a = (ofi, . . . , a*), | a| = aH h (fy , e 8a = ef (8 • • • ® e"* and 

[ fl *( e )] a = fl *( ei ) a '...a*(^) a *. 

For example, the identification between r s (Cei) <8>r j ((Cei) ± ) and T s (f¥) is done via the mapping defined 
by e\ ai <g> e V«' -> gVfai.a'), f or a ll «i G N and all a' G U^ =0 (N \ {0, 1})* . This can be iterated but re- 
member that the definition of infinite tensor products requires the additional specification of one vector per 
component which is hopefully rather canonical for Fock spaces endowed with a vacuum vector (see ||34l ) . 
Below is a notation convenient to the definition of tensor states and which allows some extensions. Consider 
the linear isometry Cj on — F s (3f) defined by its action on the Hilbert basis (<? va ) aGLr (pf*)<t 

Cje Va = , -. \ vnl a*(ej)e va = - 1 a*(ej)e va = e v(a+ ^ , (62) 
1 \a*( ej )e^\ v 11 y/e( aj + l) K J) 

with \ lj\ = 1 and = 1 . In the tensor decomposition Y s (3?) ~ V s (Cej) £g> r s ((Cej) x ), this isometry Cj 
is nothing but the tensor product — i=a* (ej) 

L v * 

Definition 4.3 Let 3f be endowed with a Hilbert basis (ej)j e m*, for j G N*, and take the isometries 
(Cj) ; e N* defined in Jt? by i62\ . For j G N*, the operator Ej is defined on Jz? (3ti?) by 

Ejp=CjpC*, VpGJ^ 1 ^). 
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For X = (Ay)y 6 flp* € ^ ([0, +°°)) swc/z f/zaf LJ=i ^7 = 1> ^ notation X.E means 

X .E = ^j^j • 

.7=1 

The operators £y and A .£ transform normal states on 1 ^ into normal states on V* 3? an d they all 
commute. After taking q>\=e\ and (fa = ei the tensor state on Y s {2f) identified with p 1 ©p 2 and studied 
above with iFi = Cei and ^2 = (Cei) x is nothing but 

E {Nt M |fl) (fl| = fif »£f |fl) (fl| = £f Ef 1 |fl) 

(A.£) JV = £ ^A"E a . (63) 



-1 ^2 - ^2 ^1 

Moreover the multinomial formula holds 

Nl 
. a 



\a\=N 

We use these notion to formulate the propagation of nontrivial Wigner measures. The Hamiltonian is 



Wick 



with (A, @{A)) self-adjoint and Qj = Q*j e Jz?(V ; 3f) ■ It is associated with the mean field Hamiltonian 

h(z,z) = (z,Az)+Y t Qj(z) 

.7=2 

and the flow (F r ) re jj in the phase space 3f . 

Proposition 4.4 Let 3? be endowed with an orthonormal basis (e ; ) ;e pj* and let the family (EA i e ^* be as 
in Definition \4.3\ Once p £ (0) is fixed pe{t) is defined by peif) = e~'e He p £ (0)e l i He . 

1) For k 6 N* and (Vi , . . . , V^) G [0, 1]* fixed such that YJg=i v £ = 1» assume that N( equals the integer part 

[veN]for£e {1, ...,£}. Then the family of states {p £ {t)) £=l / N given by p £ {0) = E^ N[ w *)|£l)(£l| admits a 

unique Wigner measure 

li, = (F,)*Mo = (FtMSj^ x ' ' ' x $%e k ) ■ 
The reduced density matrices Jg (f ) converge in Jzf 1 ( \J P 3f) to 

Jj>) 

by setting zt = F,z . 

2) Let A = (A 7 ) /G n* G I 1 ([0, +°°)) fee such that Y,°j=i Ay = 1. 77zen the family of states (pe(0)e=l/7V £' ven 
/3y p £ = (A.£') A '|i2) (£2 1 satisfies the same properties as above with 

j=l \/ k i e i 

Proof. Actually it suffices to identify the measure /Xo and to check the assumptions of Theorem ll.il at time 
r = 0. 

1) It is a simple application of Lemma |4~2l with the decomposition 



t\t) = Jjzf p )(z? p \dMz) (64) 



r s {&) ~ r,(Oi) ® • • • ®r s ((Ce M )) ®r,((c«i © • • • © Cen) 1 ) . 

In this decomposition E^ Nl >~> N k) |£2) (£2| is nothing but a tensor product of Hermite states. \ef e ){ef 1 | and 
the result is a simple tensorization of the result for Hermite states with e = ^ . 

2) The state p e (0) = (A.^jO) (Q| belongs to if 1 (V* ^) . It is therefore localized in the ball with radius 
1 . According to Proposition ^. 151 its Wigner measures are completely determined if we know the limits of 



Tr p £ (Q)b 



Wick 
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for all the b € &°°,J3f) . Due to Pythagorean summation, the measure u.q — x°° = , S 5 '^-- is supported in 
the ball of radius 1 . The estimates 



Tr 



and 



Pe (0)(b-b') wick 
(b(z)-b'(z))dii (z) 



Tr 



p £ (o)x(N)(b~br k:k xm < c p , q \b-b'y M 



3T 



{b{z)-b'{z))x 2 {\z\ 2 )d^{z) 



■CvJb-b'\ 



with the first one deduced from the number estimate (TlOb in Proposition ^. 31 hold for all b,b ! g ^p q (3f), 
p,q € N as soon as % € ^o°([0)+°°)) is chosen such that J = 1 on [0,1] . Hence it suffices to prove 
rim e ^oTr [p E (0)b wick ] = J^b(z) djXo(z) for a total set of (3f) . With the compact kernel condition, 
any b 6 Jz?°°(V / ' -2°, V 9 can t> e approximated by a linear combination of rank one operators of the form 

™ e ^l = \fM^ y \r\^^ P ^ W = P> M = «• With 

^P))™* = [a*(e)Y[a(e)f 

and p e (0) = £ A«^| e v «)^ Va |, 

|a|=iV 



we can compute directly 



Tr 



p e (0)(<z%e^)< e ^z^)) Wie * - E ^A a <fl(«) r « Vo ,a(«)V a ) 



|a|=JV 



a! 



Actually 



with a similar identity for y yields 



else, 



Tr 



p e (0)(<z®V® r )<^,Z® p }) mc * 



— <5fi v e 



With e = 1 /AT and (p,q) fixed, we obtain 



Tr 



Pe{0){{z®\e^)( e ®P ,z® p )) wick 



5^ !y e p AT(iV-l)...(^-p + l)A^. 
5^ = j^,e^){e^,z® p ) diM>{ 



□ 



We conclude with two remarks: 



• The tensorized Hermite state fi^ 1 • - JV < ? -> |Q) (£2| with N e = [XiN] and £J =1 A/ = 1 can be studied and 
behaves asymptotically like (X.E) N \Q.) (£2| . 

• When those tensor states are not Hermite states, the reduced density matrices satisfy no closed equa- 
tion and all the hierarchy has to be considered. In the example leading to d5oT l for Hermite states the 
general equation for y^ 1 ' (f ) writes 



id t7 ^(x,y) = [-A, r ^}(x,y) + J^V(x-x% Z) (x',x,x',y)-y^ ) (x',x,x',y)V(y-x') dx 1 , 

(2) (3) 

and the equation for y, involves yj and so on. . . The propagation of Wigner measures gathers all 
the asymptotic information in this case. Geometrically it is interesting to notice that if the initial 
Wigner measure is 8 s /— x 8 s j— , with Aj + A2 = 1, it is supported by a 2-dimensional torus. After 



(1)1 



.'w( 2 )rfV 



the action of the continuous flow, the support of p r remains topologically a 2-dimensional torus but in 
general deformed in the infinite dimensional phase space with no exact finite dimensional reduction. 
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4.4 Condition (PI) for Gibbs states 

For a(e) G Jzf 1 (5°), which is a non negative strict contraction: 

a(e) = 1 £a i (e)\e i (e))(e i (e)\ , 0<oj(e)<l, £ a ( -(e) < +«., 

where (ej(e)) i€ pj* is a Hilbert basis of J°, the operator r(c(e)) belongs to _S? 1 (Jf 7 ) . It equals r(a(e)) 
I^r=o i ^n( (7 ( e )) " i ^n anc ^ *h e tensor decomposition gives 

Tr[r(ff(e)) ^ni-^i) eR+ - 

Hence we can consider the quasi-free state 



Tr[r(a( £ ))] 



r(o(e)). 



It is more convenient to write 



o t (e) = with v ( -(e)e[0,+oo), 



Vi(e) + e 

and the condition Y4L1 Oi( e ) < +°° i s equivalent to Y,j=i v i( £ ) < +°° ■ 

Lemma 4.5 For a(e) = YT=i \$p^\ e i( e ))( e i( e )\ e & l (&), the quasi-free state p £ = Tr[r( ^ (£))] r(g(e)) 
satisfies 



VfeeN, sup Tr 

£6(0,6) 



if and only if there exists C > rac/z f/zaf v,(e) < C. In such a case, the quantity Tr [psN^] , k € N, is 
equivalent to 

k\ £ v(e) a 

l«l=* 

a* £ 0, w/f/z f/ze MSMaZ multi-index convention, v{e) a = EHLi v *( e ) a * • 
Proof. Consider forx e [— c,c], c > 0, the quantity 



IE 



Tr 



p £ (l + ex)<: 



1 1 jgg 

v,+e 



lli=l £ 



(e)x 



When Tr [peN^] is uniformly bounded w.r.t e <= (0, e), for all k e N it is a < ^ , °° function around x = with 



p £ (l + ex)" 



\x=Q 



■■ Tr [p £ N(N - e) . . . (N - (k- l)e)] - Tr 



P,N' 



as e -j> . 



But the first derivative is nothing but 



^Tr 



p £ (l+£x)- 



\x=0 



lv,( £ ) ; 



i=l 



which says that the uniform bound v i( £ ) < C is a necessary condition. 

Reciprocally when YT=i v i( £ ) < C, then the function l\T=i (1 — v /( e ) x ) 1 is analytic with respect to x in a 
disc of radius Rc and equals 



i=i 

which yields the result. 



co 00 00 00 

^ v^e)*)- 1 = EKE v,-(e)V) = 

;=1 ;=0 ' ' 



jt=0 



E v( £ )« 

.i«i=* 



□ 
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A Gibbs state is a quasi-free state with cr(e) = e £L ^ where L(e) is a strictly positive operator assumed 
here with a discrete spectrum: 



L(e) = £ 4(e) 4(e)<4+i(e), 



(65) 



i=1 



where the basis (ej) ; e f$* is assumed independent of £ 6 (0, e) for the sake of simplicity . There is a simple 
traduction of the assumptions of Theorem ll.il the non obvious one being the condition (PI) hidden in the 
assumption (f2|. 

Proposition 4.6 The Gibbs state p £ — Tl ^ e - E L(e)^ r(e~ £L ( £ )) with L(e) given in (|65T > satisfies the assump- 
tions ofTheorem \l.l\ ifand only if 



For all i € N* the limit lim £ ^ 4(e) = 4(0) exists in (0, +°°] . 



• If J eN*U {°°} denotes the largest element in N* U {°°} such that 4(0) < +°° for all i < 7, f/ze fwo 
conditions are verified 



J 



1 



£4(0) 



ant/ lim V — 



ee 



:0. 



Proof. First of all, writing a(e) = e~ £ (£ > allows to apply Lemma |431 with v,(e) = <v 

p -eli{e) 

Vi(e) > 



e £ ^<( £ ) > 1 — e4(e) we deduce 



l_ e -e<i(<0 



(66) 
(67) 

From 



4(e) 



Hence the uniform boundedness of Tr [p^N*] for ieN, which is equivalent to Y*7=\ v i( £ ) C implies 

inf £j(e) = K>0. (68) 

;GN*,£G(0,£) 

We now use the assumption that the family (Pe)ee(o,e) admits a unique Wigner measure po . As a quasi-free 
state, p £ is given by its characteristic function (see for example ifTSl and Q for the £-dependent version) 

Tr[p £ W(/)]=e 

But the Wigner measure is characterized by its characteristic function 

G(|) = / e- 2iKS{ - z ® dno{z) = limTr [p e W(V2^) 
By taking | = <?,-, ieFf, this implies that the limit 

lime 

£^0 

exists in K. With the constraint d68l l there are two possibilities: either lim e _ > o4(e) = 4(0) € [jc,+oo) and 

G(e,) = e l: i {0 ) or lim £ ^o4(e) = +°° and G(e;) = 1 . After recalling that the 4(e) are ordered and by 
introducing the index J like in our statement, we get for | = Y^Li € .2° 

.2d Jill 
G(|)=e Li=1 Wl . 

The measure po has to be the gaussian measure 



J 

p = x 
(=1 



4(0) -e,(o)kif 



L(d Zi ) 
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Our assumptions imply that the integral \z\ 2 djUo(z) equals 



z\ 2 d^l {z) = limTr[p £ N] 

«J° £-5-0 



After Lemma 1431 we know that 



7 1 
I 



i(0) 



lim Y Vi(e) = lim y 



-eti(e) 



which enforces the two conditions d66l i and ( |67] > . 

Conversely assume that all the conditions are satisfied. Reconsidering the final argument in the proof of 
Lemma POI savs that the function 

ft (l-Vi(e)x)- 1 

i=J+\ 

converges to 1 in a given neighborhood of x = . Hence 

/ ee -eti{e) \ " 



lim Tr 

£^0 



p £ N s 



limit! 

£^0 



I 

\a\=k, 
a,=0 for i>J 



1 



k\ 



I 

a,—0 for i>J 



£(oy 



which is easily checked to be equal to Jj» \z\ 2k dno(z) . □ 
In the Bose-Einstein condensation of the free Bose gas in dimension 3, considered in [7|, the first 
eigenvalue is tuned so that ^i(O) <E (0,+°°) and all the other eigenvalues are such that £/(0) = +°° . The 
condition which fails and gives rise to a physical example of dimensional defect of compactness is d67l i . 

4.5 The Hartree-von Neumann limit 

Let po be a non-negative trace class operator on L 2 (R d ) satisfying Tr[po] = 1 and let 

p® A, = p ( g ) ...®p. 

Consider the time-dependent von Neumann equation for a system of particles 

f id t p N (t) = [M N ,p N (t)) 



Pn(0) 



Po 



(69) 



with p#(f) is a trace class operator on L 2 (M </ ) I?)A ' ~ L 2 (M. ). Here is the Hamiltonian of the particles 
system 



N 
i'=l 



KJ 



with A is a self-adjoint operator and V € L°°(R ) real-valued satisfying V(x) = V(—x). As will appear in 
the proof, more general interactions could be considered in the spirit of Theorem ll.il but we prefer to stick 
to the usual presentation for an example. 

The next result concerns the limit of the von Neumann dynamics d69l in the mean field regime — > °° 
already studied in [ 10] [9 1. We shall see that although the particles are not assumed to be bosons, our 
bosonic mean field result apply to this case due to the symmetry of the tensorized initial state pQ N . 

Proposition 4.7 Let (p^(t)) denote the solution to ( 1691 ), and consider the trace class operator ajp(t) G 
££ 1 (L 2 (R kd )) defined by relation 

VBGJ£f(L 2 (R M )), Tr\a { N k \t)B 

Then the convergence 



Tr 



PN (t) {B (8> I L 2 Qg,d(N-k) ) ) 



lim a, 



(k) 



N- 



N 



■■Pit? 



(70) 
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holds in 3? l (L 2 



or alltdzW and when p (t) solves the Hartree-von Neumann equation 



id t p(t) = [A + (V*« p(r) ),p(f)] 
P(0) = Po, 



(71) 



with n p (x,t) :— p(x;x,t) . 



Proof. The proof will be done in three steps: Bosonisation, Liouvillian and mean field limit. 
Bosonization: The phase space that we will consider is not the one particle space L 2 (W l ) but 

iT = ^ 2 (L 2 (R l1 )), 

the space of Hilbert-Schmidt operators on L 2 (R d ) . It is endowed with the inner product 

(ffli,a)2>^ = Tr i 2[a) 1 *C()2] 

where Tr L 2 [.] here denotes the trace on L 2 (K ) and ©jf is the adjoint of a>\ . 
The cyclicity of the trace leads to 



Tr 



(L2)« 



p N {t){B®I Ll{wl(N - k)) ) = (^ N (t),(B(x)I L2{Kd{N _ k)) )^ N (t))^ m 



(72) 



with m N (t) = e - m ^ y /p^' N e m « . 

The important point is that at time t = 0, IV (0) = y/po® , is a Hermite state in 3f and that the evolution 
preserves this symmetry so that 



With any bounded operator B : L 2 
by 



> L 2 (W lk ), the action by left (resp. right) multiplication is defined 



L B (resp. : \J k 2f -> \j k & 

(0® k i — ^ y k {B(0® k ), (resp. ^((B«*2?)) , 

where is the orthogonal projection from Cg^iF onto ^ ■ Since (©^oej" is a total family in 
this defines a bounded operator L B G j£?(V* ^) (resp. fl B G .if (V* ^0) such that L* B — Lb* (resp. R* B =R B *) . 
When ,x k ,y\ , . . . is the Schwartz kernel of B G j£f (L 2 (R d *)) , Lb (resp. R B ) is the left (resp. 

right) multiplication by the operator with kernel 



k\ 



y 52 ^( Jf ff(i))---»-fff(t)))'ff(i)i---))'ff(jfc))- 



Hence the trace d72l equals 



Tr, 



(L2)S 



PN{t){B^Ij2.^d(N-k))) - (^Nit) ,L^ mie>{N - k) ^¥ N (t))y N ^ . 



With an operator B G -2?(L (K *)), we can now associate a symbol 



b B ((0) = ((Q® k ,L B (0® K )y ks , = Tr (L2)S * 
Since £r fi g>/®(Ar-*)i is nothing but L B \J IyN- k ^ we get 

{N — k)\ 



Tr 



(L2)S 



PiY {t)(B® I L 2 ( R rf(JV-i) ) ) 



Me* 



V 2? ' 



1 



Liouvillian: Let us now determine the appropriate Hamiltonian H £ of this problem which is actually a 
Liouvillian. The map 



3t^te i,A coe itA 
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defines a continuous unitary group on 3f with a self-adjoint generator 

Z A : 2T -> 

0) i ^ [A,o]. 

The interaction is a bounded self-adjoint operator Q:\f 2 ' V 2 3f given by Q = \ (L v -R v ) E Jz?( V 2 3?) 
and we associate the symbol Q(co) = (<B 02 , go® 2 ) . For any CO E ^ the kernel of Qco® 2 e\/ 2 3? is given 
by 



(2co® 2 )(jci,}'i;x2,y2) = 2 V ( x i -*2)»(*i,yi)<0(*2,;y2) - ^Cvi -y2)(o(xi,yi)a(x2,yz) ■ 
After introducing the Hamiltonian 

H £ =dT{Z A ) + Q mck , 
acting as a self-adjoint operator on T s ( 3f), we get for © e V £^(dr(£ A )) , 

e l H £ ®= [IHW,©] with e=l/N. 

This implies 

Mean field limit: The initial data p e (0) = l-y^A) ) (y'Po 8 ^! is a Hermite state which fulfills the assump- 
tions of TheoremQTT] with 

^ = 5 % ■ 

The classical energy associated with the Hamiltonian H £ is 

h((o) = ((o,£ A( o)^ + ^((o® 2 ,(Lv-Rv)(o® 2 }^ 

and the mean field flow F f is nothing but the one given by 

idfCO = d a h(a>) = \A,co] + (V*n 1 (0 ) to - co (V * n 2 m ) , 

where V *n' a are multiplication operators and n l m {x) = f K d \co(x,y)\ 2 dy, n^iy) — J K d \(o(x,y)\ 2 dx when 
G)(x,y) denotes the kernel of CO . Beside the invariance |F t (fl))| %' = \<&\x an d F,(e~' e co) = e~' e F,(co) , 
the flow F f also satisfies 

F f (fl>*)=F,(fl>)*. (73) 

Thus previous equation becomes equivalent to the Hartree-von Neumann equation (T7T~b with p(f) = Co(t) 2 
when co(0) = -^Po . The Theorem ll.il says 



Vbe0> Kk {2f), limTr^ 



\V N (t))(V N (t)\b 



Wick 



In particular when B € S£ [L (R )), this implies 



lim Tr 



p W (f)(B(g)/ L 2( R rf(Ar-t) 



Tr 



L 2 (K*) 



p(t)® k B 



This proves the weak convergence in d70l >, but since it is concerned with non negative trace class operator 



and Tr 



a N ' (0 = 1 = [p (f ) 8A: ] the convergence holds in the Jz? '-norm. 



We end with three remarks: 

• When p is a pure state, the result of Proposition l4.7l is the same as 



□ 



When p is not a pure state the Subsection 14.31 has already shown that one has to be very care- 
ful with tensor products. Actually p®" £ Jz? (gr JT) commutes with the symmetrization projec- 
tion J^y (or the antisymmetrization &/n for fermions) but the corresponding states in Jzf^V 3f) 
(resp.JSf 1 (A JV 30)are 

.y N p m y N (resp. *f N p® N £/ N ) . 

But as shows the formula Tr^^p)] = Ulea(p) izj ( res P- Tr [ r «(P)] = Ulea(p) j^), the trace 
of S^isip® 1 * 'S*n (resp. ^p® N \s/n) converges to as N — > °° . We leave for subsequent works, the 
question whether normalizing these states would lead to the same asymptotics as in Proposition l4.7l 
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• We recall that a tensorization based on the tensor decomposition of Fock spaces in Subsection |4.3| led 
to the evolution of Wigner measures which cannot be translated in terms of Hartree-von Neumann 
equations. 

Acknowledgements: This work was finished while the second author had a CNRS-sabbatical semester in 
Ecole Polytechnique. 
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